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ABSTRACT 


This thesis is concerned with the development of a non- 
Pineda Dinary-distillation model, to include the hydraulic 
effects of variable volume tray holdup. Verification of this 
model is made by comparison to open loop data obtained from a 
pilot plant column. In addition, other limiting assumptions 
are consecutively applied to the model, and effects on model 


behavior observed. 


A literature survey presents various types of distilla- 
tion models and their applications. A list of common model- 
ling assumptions is given along with effects upon model com- 
plexity and order. A description of workers' treatment of 


Special stages such as reboiler and condenser is included. 


Mass, component, energy balances, and a linear hold- 
up correlation were used as a basis to develop a nonlinear 
model of a ten stage binary column. Two metheds of numerical 
solution of the resulting forty equations were developed. 


Simulations were made for reflux and steam flow disturbances. 


The model was simplified to include the assumption of 
constant tray holdup. Further limiting assumptions were 
consecutively applied and simulations were run for steam, 


reflux and feed flow disturbances. 
Verification of the constant tray holdup model was 
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achieved by comparison to open loop transient data obtained 
from a pilot plant column. Very good agreement was shown 
with heat loss and efficiency used as fitting parameters. 
Comparison of the constant tray holdup model to the variable 
tray holdup model and other more simple models indicate the 
effects of various assumptions upon model behavior. Simula- 
tion results from the variable tray holdup model with three 
different holdup correlations are shown. It is shown that 
the choice,of a holdup correlation can significantly affect 


the simulated column dynamics. 


The final model formulation was found to lend itself 
easily to modification, and handled the variable holdup con- 
sideration well, providing a stable solution for a variety 


of cases. 
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INTRODUCTION 


It has been stated by prominent chemical engineers and 
designers that process dynamics and control should be con- 
sidered in the mmrtial design of a process., Such a design 
‘philosophy would include the development of a mathematical 
process model to aid in improving process operation on the 
basis of dynamic and control considerations. This approach 
has not been widely applied, with the control system gene- 
rally being designed after the equipment design has been 


finalized. 


The desired design philosophy.must combine the mathe- 
matical representation of the process system with tentative 
control schemes. Simulations of the controlled system be- 
havior would lead to improved process designs. The applica- 
tion of this approach has been limited due to the difficulty 
of understanding the process physically and developing ade- 
quate mathematical models at reasonable cost. It is for 
these reasons that such a design philosophy has been success- 


fully applied only to the more simple process systems. 


Distillation, however, is an exception, being a rela- 
tively complex process. It has been the subject of many 
studies in regard to dynamic modelling, as it 1s one of the 
oldest and most widely used process units in industry. Models 
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of varying complexity have been derived, some having been 
VeriLled by Comparrson to actual pPlot plant dataw” iThese 
models include limiting assumptions that determine their 
rigor or completeness. A complete description of simplifying 
assumptions commonly used in distillation modelling is pre- 


sented in a later section. 


A major assumption which is made in deriving any distil- 
lation model concerns holdup and hydrodynamic lag. The pri- 
mary purpose of this work was to eliminate any assumption 
concerning holdup or hydrodynamic lag, augmenting a fully 
non-linear model already incorporating a minimal number of 
assumptions. This proposed model describes open loop tran- 
Stents of tray and terminal compositions for a variety of 
disturbances, including changes in feed, steam and reflux 
flows. The model is verified by comparison to open loop data 
Crom aessinch diameter, LO stage pilot plant distiltation 


column using a methanol-water feed. 


Models developed previously have given rise to favorable 
prediction when compared to experimental data, or have been 
used to simulate successfully a column under a new control 
scheme. However, few comparisons have been made between 
models of varying complexity to determine the error introduced 
by particular assumptions, reduction of order and/or lineari- 
zation. Thus, the objective of this work was to use a com- 


pletely non-linear experimentally verified model and then 
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modify the model by more restrictive simplifying assumptions, 


determining their effects upon simulated column behavior. 
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CHAPTER 2 


LITERATURE REVIEW 


2.1 ASSUMPTIONS 


Zabaleaulntroduction 


There have been two distinct approaches to modelling 
distillation and other processes. The first and most elemen- 
tary is fitting a transfer function model of predetermined 
order and complexity to open loop responses between desired 
input and output variables. The second, and more flexible 
method is the application of known physical relationships such 
as mixing equations or mass and energy balances to sections 
of the process and predicting open loop responses. The com- 
plexity of the models that would result from such an approach 
is determined by the assumptions employed. The 
most common assumptions made in distillation modelling are 


listed below: 


1. Negligible vapor holdup and dynamics - changes in 
vapor rate are assumed instantaneous, and vapor holdup is 
considered negligible relative to liquid holdup. Use of such 
assumptions, which are normally justified, significantly re- 
duces the number of equations to be solved. 

2;:-Molar,-mass-or.volume holdup constant + the specific 
assumption chosen depends upon the extent of the simplifica- 


tion desired. This simplification has been used in almost 
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all distillation models to date. 

3. Liquid dynamics negligible - with this assumption it 
@esnconsidered that’a’ disturbance in liquid flow rate will 
propagate through the column instantaneously. This may be a 
good assumption if the composition anpnieines are very slow, or 
if there is only a slight dependence of holdup upon liquid 
flow rate. 

4. Constant pressure operation - simplifies calculations 
in that the effect of changes in pressure on equilibrium or 
enthalpy considerations need not be considered. Tower pres- 
Sure is under control in most processes. 

5. Well mixed liquid and vapor phase in each stage - 
only stage efficiencies are required to establish liquid and 
vapor compositions on each tray. 

6. Adiabatic operation or other heat loss assumption - 
these assumptions are only applicable if energy balances are 
used in the model. Adiabatic operation may be valid for well 
insulated or enclosed columns. Other common assumptions are 
those of constant heat loss or heat loss varying with compo- 
Sition and/or temperature. 

7. Constant efficiency or some function of operating 
conditions - no simple correlation for predicting the effi- 
ciency of small columns exists so efficiency is used as a 
parameter to fit experimental data at steady state and/or 


during transient operation. 


: ; nila 
. + : 
- i &, . 
| t a 
ened / 
a Da , . 
) ; 4) lies = 7 7 oa : 
ay: . — 4 aan 7 i i : 
$75b 0 ae aiebem sorts 
es cdiw ~ eidey Seam =e thalia 
“ fe gd 
i aot? tive sane ai & tsi i 
7 << ‘an “n met he 
a) ; (Leuns tet isteaL ani any aipeiat ’ 
"a 
: ie ; a) 1] 7} yi 
2 - ve Page » e « 
; 
* { 
i 2 wc 
- L= pe rc im AB 4c os oe S23 
( 7 
$4, ey ie «J wie = eb 
bawWe DAGSe «Sey 
4 | 
: Wi 
j 2 af — wT) no = ie it =e 5; Caius Li ia’) * ‘G 
SL : = a bpres iG "ee tine aan 
vE = = na hae k 
& Me f f mae 
= if. Wirs pas at 2. Jar i ha 0 ae.) S aie = ei. a 
af siehelse wpasce SS Bluest tags me NO .2R4 sro] 
ay —— 
iiow 2s 57 ew Se) Gaia teas SESE CSEID. 
1 7 i 
S46 G2 Ioticas patwis SSte oeiaiics Suet 


—ietys orl phi sev Ribas 
a a ee. 989 
; | 


: we 2 
a a. 
oan On 

; : 
: : > ¥ i, 
Tr¢@: oon a 
oe. ne 


12 4zeel ’ 24 Sri + SESS 
~>;' 


FOS 3 — > ~ 


7 eaeh 


¥ 


a 


Ada 
‘ob Lig 


set 


8. Constant molal overflow - the liquid and vapour flow 
rates are constant in each section of the column. This is a 
severe assumption and usually only introduced with minimum 
prior knowledge of column operating conditions. 

9. Perturbation theory applies - only small excursions 
from steady state will be experienced which allows equilibrium 
and enthalpy relationships and/or conservation equations to 
be linearized. 

10. Heat capacity of tray metal and/or column walls 
negligible - energy storage considered to be small compared 


to magnitude of other terms in enthalpy balances. 


As shall be seen in later sections, all modelling work to date 
has included the use of at least two of the above assumptions. 
The accuracy of models predictions using certain of the above 


assumptions will also be considered in later sections. 


2.1.2 Holdup and Hydrodynamic Lag 


To a great extent in the literature the effect of varying 
holdup and hydrodynamic lag have been considered to be negli- 
Sible. “Physically, "trquid flow from a’tray is determined by 
the head of liquid on the tray, the amount of frothing, and 
to a certain extent, the pressure differential and vapour flow 
to the tray. The exact nature of this relationship is unkown. 
For conservation of mass, an increase in input liquid flow to 


a tray would result in an increase in liquid outflow and 
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possibly an increase in liquid head on the tray. As in other 
reservoirs with a restricted outflow, the relationship may be 
a lag equation 

L(s) ak 


LAS) ; tse (2.1) 


and the holdup dynamics would also be in the form of a lag 


NSD) fea Nok 
L. (s) Ti rSel (Ain) 


Proposed models reported by Stainthorp and Searson’, Scho- 


3 


lander’ and Waller” have used expression (2.1) to relate 


liquid flow rates. 


Many linear models assume constant liquid flow rates 
and holdups (L(t) = Lo: WT(t) = WT) upon which composition 
transfer functions or matrices are based. Such is the case 
with models developed by Osborne et at ae by Sittel tand 


Fisher? and by Rosenbrock®’ ’, 


Holdup and hydrodynamic lag are treated in a different 
manner in the more complex nonlinear models. Liquid and 
vapor flows at any time are determined by mass and energy 
balances on each tray, taking into account accumulation or 
holdup changes due to density or composition changes. This 
eliminates liquid dynamics altogether with flows changing 
only with composition effects or heat effects. This approach 


9 
was adopted by Huckaba et al” and by Rafal “et al. 
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in their studies. 


As stated in section 2.1.1, the most common assumption 
made with respect to holdup is to assume the holdup to be 
constant on amass, molar or volume basis. Dee ean ae has 
altered the model of Huckaba et Sang to correct tray holdup 
for changes in liquid and vapor rate. Holdup was updated 
each integration interval, but changes during each time step 


were not considered when solving for new liquid and vapor 


flow rates, resulting in neglect of liquid dynamics. 


esc low OL, Liquim. room one, tray to-another ius through 
a downcomer, its length approximately that of the tray spacing. 
One cits Legion, liquid does not come into contact with the 
vapor, and should be modelled as a pure time delay in liquid 
flow. This approach was used in the models of Peiser and 


eroverus and that of eon 


fos MOllObv ous Lf there sa basis that: should be vem-— 
ployed regarding the use of particular assumptions since good 
results were obtained by previous workers using a particular 
model for a specific case. Only upon comparison of experi- 
mental data with simulation results from different models can 
a judgement be made. Even the prediction of actual holdup 
values for bubble cap or sieve trays involves considerable 


uncertainty. 
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2.2 TRANSFER FUNCTION MODELS 


The earliest distillation column models tended to be 
the simplest, and transfer function models fall into this 
category. Transfer function models are easily obtained, by 
fitting to open loop data or by simple mathematical analyses 
of the system. Transfer function models only characterize 
nonlinear plants adequately provided there is no strong non- 
linear system behavior. The distillation column is a non- 
linear system; however, for small excursions from steady 
state, the transfer function model may be representative. 
However, with the linear model and a highly nonlinear system, 
modelling error became important, and with the wide availa- 
bility of digital computers rather than analog computers, 


other system representatives were derived. 


A process transfer function relates the dynamic res- 


ponse of an output variable to a disturbance of the input 


variable. For a process transfer function is given below. 
at C(s) (203) 
G(s) = a(S) 


The time domain response is obtained from the inverse 


Laplace transform according to 


Ci mers tete) us) 4 (2.4) 
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In a distillation column, the variables of interest or 
output variables are usually the terminal compositions or 
compositions of the product streams. The disturbance variables 
are usually the feed flow, feed composition and feed enthalpy. 
The relations between heat input (steam flow) or reflux and 
bottom composition or top composition are also important for 


Gontrol studies. 


Wahl and Harriot?? attempted to predict the time con- 
Stants of these transfer functions for a 15 plate methanol- 
water column by fitting open loop data and by using previously 
developed correlations based upon a recirculation rate anda 
characteristic time constant (column parameters). The cor- 
relations predicted both zeros and poles from steady state 


Sonds-crvon's:; 


Transfer function models relating certain tray tempera- 
tures to reflux and steam flows, determined empirically for 
+10% step changes in flow rate, were used by Vinante and 
Luyben!? EOGedecoupsena control fot ja, dist: Llation column. 
Luyben and Tyreus-®, using a first order hydraulic dynamics 
transfer function model fitted to experimental data, evaluated 
sidestream composition control. Other workers who have derived 
transter function. .modelsware.. Sitteland aha ee Paul 


et oat and Waller, the latter assuming that flow lags and conm- 


position lags can be treated independently. 
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Another approach has been to simplify higher order or 
more complex models by fitting their predicted responses with 
Emenster functions. Typical’ are the studies of Izawa and 
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Morinaga”™ ‘and Beaverstock and Harriot*”>. 


2.3 STATE SPACE MODELS 
State space models may be represented by Equation (2.5) 


aha Xt 0 (2945) 


moc eeR useaematraxson coefficients and’xpista state vector. 
The state vector is a convenient way of representing impor- 
tant system operating variables, individually or as linear 
eombinations of sother variables.! ‘Thesmatrix A isja matrix 

of coefficients, predicting trajectories of the system through 
state space. The other terms on the right hand side of Equa- 
eLOneW2s5) aremcomprisedisimitarlyeot: aymatrixsot coefficients 


and a vector which may be disturbance or control terms. 


State space models may be derived by mass and/or energy 
balance considerations. Each tray in a distillation column 
has composition, energy and total balances associated with 
it, as well as equilibrium relationships and holdup conside- 
rations, which may all be used in expressing a state space 
model. These relationships, however, must be linear or 
linearized. The state vector may be plate compositions and/ 


or enthalpies, augmented by liquid or vapor rates or holdups, 


ie 


depending upon model rigour or complexity. 


Scholander“ used a model based on material balances to 
calculate feedforward and feedback matrices for bottoms 
composition control. Rafal and Stevens” derived optimal con- 
trols based on quadratic performance indices for a distilla- 
tion system using a state space model. A model by Clear*? 
used K values to predict equilibrium conditions, and heat and 
mass balances to describe column dynamics. Binder and 
Calvillot® used a state space model to test a tracking con- 
trol scheme and an adaptive-predictive control scheme. The 
derivation of these various control laws, and the use of 
other standardized techniques for model reduction are all 
matrix operations, and are possible because of the state 
space models linearity property, and the form of the model. 
Terris’ used a matrix technique to reduce a state space model 
to order 2 from 121, comparing time constants and gains of 
the reduced model to industrial data. Simplification of this 


kind may reduce computation or complexity in application 


but with some loss of accuracy. 


Pike and Thomas?3 obtained a fourth order state space 
model by discretizing Laplace transfer functions and conver- 
ting to a matrix form. With mathematical techniques appli- 
cable to matrices but not to Laplace functions, they calcu- 


lated an optimal control strategy with two variable control. 
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the tinedricy and practicality of the state’ space models 
allow a variety of control strategies to be rigorously cal- 
culated and implemented. Unfortunately, such a model has 
found limited application since a distillation column is a 
nonlinear system. Further, without fitting of some parameters 
such as efficiencies and/or heat loss, even steady state 
agreement between calculated and experimental conditions is 


unlikely. 


2.4 NONLINEAR MODELS 


A nonlinear model consists of a system of nonlinear 
differential and algebraic equations. In practice, discrete 
methods for the solution of nonlinear differential equations 
are as applied to linear equations with changing coefficients; 
the solution for the nonlinear distillation system is ob- 
tained by solving the composition balance equations, of the 
£Orm given by Equation’ (2.5), and updating the A matrix for 


Changes in the state vector x. 


The nature of the updating technique is again dependent 
upon the nature of the simplifying assumptions chosen. When 
the nonlinear model is derived from mass, energy and compo- 
nent balances, such as in the work of Huckaba et bgouty the 
system of differential equations involves the liquid and 
vapor flow rates, vapor compositions and stage holdups. These 
may be established via a number of routes, and using various 


assumptions such as constant molal overflow and 100 percent 
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tray efficiency as in the work of Osborne et a ibe Distefano*4 
solved for the composition coefficients, as was done by 
Huckaba and coworkers, but included the effect of a dead time 
in. liquid flow from tray to tray and the effect of changing 
volume holdup on the tray. Svrcek*® modelled a distillation 
column similarly but assumed constant volume but variable 

mass of liquid holdup on the tray. Peiser and Grover*t used 
mass, energy and component balances to derive a model based 

On a 4 component system. The numerical techniques used to 


solve the system of equations differ from case to case. 


The applications of non-linear distillation modelling 
have been to test proposed control schemes, or have involved 
comparison of dynamic simulation to experimental data. Few 
comparisons have been made between linear and nonlinear 
models, or between nonlinear models based on different as- 
sumptions. One such example can be found in the work of 
Rafal et alscs who compared two methods of linearizing phy- 
Bical property relationships? but found. no Significant: dift- 


ference. 


To compare favorably with experimental data, a certain 
amount of fitting has been required. Reliable correlations 
for holdups or efficiencies are lacking, forcing workers to 
adopt a fitting approach. Distefano*4 developed his own 
empirical holdup correlation for lack of a better one. 


2 . . . 
Beaverstock and Harriot’> reported in their pilot plant 
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modelling effort that the Francis Wier formula predicted 
holdup values larger than measured by a factor of 2 and hy- 
draulic time constants one third the measured values. Other 
workers assumed 100 percent efficiency when methods of effi- 
ciency prediction were available but supposedly did not 
epply vin their particular application. The weakness of deve- 
loping more complex models has been a lack of consistency in 
predicting the values of the increased numbers of parameters 
required in higher order models. 


( 


2.5 SPECIAL STAGES 


On a macroscopic basis, all stages may be treated equally, 
as the general equations describing stages include special 
berms such as heat input or liquid and vapor drawoffs. The 
reboiler and condenser however, being terminal stages may be 
termed special as in both cases there appears a time varying 
heat term usually orders of magnitude greater than in "ordi- 
nary" stages. It is this heat term and special mixing con- 


siderations that have given these stages a special status. 


Heat input to the reboiler involves the transfer of heat 
EEOm one medium, such as steam, to. the reboiler liquid. [In 
most studies, reboiler heat duty, rather than condenser duty, 


is taken to be a specified process degree of freedom. 


Reboiler modelling involves heat input prediction and 
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some mixing assumption. Peiser and Grover?’ 


used enthalpy, 
mass and component balances about the reboiler to obtain 
transfer function models relating liquid composition in the 
reboiler (assumed perfectly mixed) to steam and liquid in flow 
disturbances. Heat input was taken as a linear function of 
average temperature differential, and was varied with time. 


10 also used heat balances to solve for heat 


Distefano et al. 
input, based on steam temperature, and reboiler liquid compo- 
Sition. As in the work of Peiser and Grover, complete mixing 
was assumed. ‘Based upon actual pilot plant measurements, 
Distefano et al. considered errors of 5 percent due to corre- 
lation inaccuracies to be reasonable limits for heat input 
mreciction. Sproul and Gerster, using composition balances 
only but including entrainment effects solved for a composi- 
tion profile at various points in a reboiler system. They 
found that dynamics were changed by varying the entrainment 
Bactor. Svrcek *° showed mathematically that for his pilot 
plant column, heat transfer dynamics were rapid compared to 
the composition dynamics, so were considered instantaneous. 
Heat input has also been considered constant, as in all trans- 
fer function models, and as in the matrix model used by 


Plear. ~~ 


The condenser is a special stage in that a large quantity 
of heat is removed. All mass entering this stage is condensed and 


withdrawn as product under level control except the liquid 
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that is returned to the column as reflux. Reflux rate and 
enthalpy are the two process degrees of freedom that are 
specified. Thus the condenser is characterized by one com- 


position differential equation. 
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CHAPTER 3 


MODEL DEVELOPMENT 


3.1 THE EQUATIONS 


ee.) Tne Distillation Column 


A distillation column may be considered to be a series 
of mixing stages, in which the vapor and liquid streams 
enter, and leave accompanied by interchange of mass and energy 
between the phases. Condensation and flashing result in par- 
tial or total equilibrium between the vapor and liquid 
phases. The degree of foaming, frothing and mixing in each 
stage determines the amount of interfacial contact which in 
turn limits the extent of equilibrium achieved. On a micro- 
scopic scale, rigorous mathematical characterization of these 
physical processes is lacking, so the approach to modelling 
taken in this work was macroscopic in nature. It 1s assumed 
that complete mixing, mass transfer and heat transfer equili- 
brium states were instantaneously achieved, and that overall 


balances were satisfied. 


Insta, cdistiliation= column, although there are many 
changing variables and properties, such as diffusivities and 
mrscosities etc...) the basis of column behavior can be’ ex- 
pressed by equations in "4n" variables, (where "n" is the 
number of stages). These are the mass circulation rates 


(both liquid and vapor flows from each stage), the mass 
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FIGURE 3.1- A Theoretical Stage 
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holdups on the stages, and the liquid compositions on the 
Stages. Other variables such as density, enthalpy and accu- 
mulation rates follow directly from these base variables. 

The governing equations (including liquid feed and sidestream 


drawoff terms) are: 


mass balance 


qwT 
QO = Acaae 7 Lisl - Li + V wi nee - S + FP (Sak) 
component balance 
qwT x, 
Oe Git Pg es aL i Bes "7 a Sal 
38 Pee nT aon Bs Pn Fn (3.2) 
enthalpy balance 
qwrt ha 
ws ate er D1 nt ie Ln i. Tee 
+ Vital - Soi a F Ben - Q (337057) 


The column is said to be at steady state when the rate 
Om Changetoc) concentrationsionreach tray, ivandathe rates of 
mass and energy accumulation, that is the left hand sides of 
Peustronceds. Wi, Woo) rand: (3.3), are identically equal to 


zero. 


When the system is not at steady state, the accumulation 
terms are not zero, but some instantaneous value, such that 


all balances are satisfied. Since these balances then 
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compromise 3n equations, a total balance, an energy balance, 
and a component balance for each tray, a further n indepen- 
dent equations will be needed to solve for the 4n base un- 
knowns. These equations are algebraic in the form of holdup 


eorrelations. 


Wlen=carL wee bNVe +c (3), 4) 
n Ig Mak lah” jel n 


The values for ane bo Co which may be functions of 
liquid and vapor density are arbitrary. No correlation 
exists to date which adequately represents experimental data 


for a wide variety of columns and operating conditions. For 


oo = 0, the equation reduces to being a wier type correlation. 
This set of equations was chosen for several reasons: 


ile The equations are expressed in terms of basic chemical 
and physical properties and flows that must be satis- 
fied at any instant of time. 

Zee There are no assumptions regarding the nature of the 
liquid dynamics, only that the holdup correlation is 
valid. 

oe The equations are general and so are applicable to any 


column or set of equilibrium stages. 


For a binary mixture, thermodynamic data or correlations 
necessary to determine the coefficients in equations (3.1) 


to (3.3) can be expressed as 
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AL = A(x) 37560) 
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Since the pressure in the distillation column is con- 
sidered constant, no allowances for pressure dependence have 
been incorporated into the thermodynamic or equilibrium re- 
Presentations. | Expression (3.7). includes not only the basic 
equilibrium data, but takes into account the efficiency of 


contact. A Murphree vapor efficiency defined by 


a ve ee 6 
EA = ital Ne a} yee 
Yn-1 Yn 


is applied to the equilibrium vapor composition to obtain the 


actual vapor composition according to 


= << 
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Bee2iaothe Reboiler 


The reboiler, unlike the trays, contains a relatively 
large liquid holdup, and has no input vapor flow but has a 
large heat input. Heat input is by a shell and tube heat 
exchanger, with condensing steam providing the energy to 


flash the methanol-water solution. 


The heat transfer depends on the overall heat transfer 
coefficient, the steam temperature and the reboiler liquid 


temperature. 
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cy ea (31:0) 


The heat transfer dynamics are composed of dynamics of 
the temperature on the methanol-water side of the reboiler, 
combined with steam side temperature dynamics. The methanol- 
water solution temperature is a direct function of solution 


composition 


Cee ny) Gea 


so temperature dynamics are identical to composition dyna- 
mics. Reboiler composition dynamics are included in the 
model as for any other stage in the distillation column, by 
calculating a component balance. The reboiler equation is 


dx, 
rs - L,x, - V (3h 12) 
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By a total mass balance, and assuming constant (approxi- 
mately) reboiler holdup by level control, the following 


equation 


indicates that the reboiler time constant or reboiler resi- 


dence time would be 


WT, 
1 = 


R L. 
This time constant may be shown to be several orders of 


magnitude larger than the time constant on the steam side of 


the reboiler by examining a steam mass balance. 


es 


96 aphmboyh Io Hesgogeo StS aot ae 


sé i. 
ig coma) lo 
al Voll = 


pati sc. ebis ct a RPA epee ete heh H aed Soups: 


B | f y - 
~ises tei ‘odS* bind gnascarmia’ 3 aleaces: 
io, 3 z 
a ‘v- ; 
eI 


a t ~ 
: ae = — ye a ee tw _< = a m7 
fonidgam’ set | .2Stmecy). <ul Serene" Sele nage ea 


' 
wal 
' 


ait Yo stae we? Be -~Tewets ae See ma ed 


7h 
i 


tm) Halts aut <cumelb 2 Sk: 4a 


~ 
‘ 
\ 
cee Peet Co poTtne.s) SS" Boe ave 
— - ' 7 ee 7 
il Trek. | ee i \ ) oF a ae —- oe a te 
foOLo ey ~ » &SreJe i 


24 


The energy balance on the steam side enables solution 
for the steam temperature, and is expressed in equation (3.13), 


neglecting steam vapor holdup terms. 


aWT.A, 
Seawea — S;Ho, — Soh. ad Q (esr. es) 
BupStitucion Of equation (3.10) (and S; = S)) gives 
dh. 
WT, at = S, (Ho; = Ag) = UA(T.-T,) 


The enthalpy of the saturated condensate may be approxi- 


mated by a linear function of temperature. 


Ag = A + BT, 


Substituting and gathering terms in T. yields 


dT. 
eS Ses oe = - UA(T_-T 
& Se S; (Hg; A BT.) UA ( S R) 


or rearranging 


> se san. ak 
BWT, aT coe S, (Ho; A) 


S R 
(S,B+UA) dt S (S.B + UA) 


(3.14) 


Expression (3.14) represents a first order system with 


thermal time constant 
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which when compared to the reboiler composition time constant 


is found to be several orders of magnitude smaller. 


Equation (3.14) could be integrated to solve for changes 
in Tos but in this work, due to the significant difference 
in the magnitude of the time constants, the steam temperature 
was assumed to reach its steady state temperature instan- 
taneously. Neglecting the transient term in equation (3.14) 
yields the following equation for steam temperature. 
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To calculate the heat input at any instant, the steam 
temperature was calculated using equation (3.15) and the re- 
boiler solution temperature determined from the reboiler com- 
position. These two temperatures were then used in equation 
io LOyealong with thevoverall ‘heat transfer costficient co 


calculate the heat input Q. 


B.a.s Ine Condenser Model 


Liquid level control in the condenser operates with a 
hydrostatic head type of measurement rather than direct 


measurement of level, so either a density change due to a 
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Change in composition of the fluid or an actual level change, 
will cause a deviation from the level controller setpoint. 
Since mass holdup and not necessarily level is being main- 
tained constant, the total mass balance will not contain un- 
steady state mass accumulation terms. On this basis, the 


mass balance may be expressed as 


i Beet) (3.16) 


The composition balance for the condenser may be written 


as 
AX, 4 
Ome AL eos, . Seal 10 
he 
mm 30 
Ue 7 Ole O-ao (2) 


Equation (3.17) was used directly in the integration 
step of the model solution, to give the top product composi- 
tion. The reflux and distillate streams are at the same 


composition, feeding from the top product reservoir. 


The amount of heat removal from the condenser was not 
included in the model, but the energy balance degrees of 
freedom were defined by specifying the enthalpies of all 
streams entering and leaving the column and by specifying 
heat input to the reboiler. The enthalpy of the reflux 


entering the column was taken to be the enthalpy of saturated 
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liquid at the top product composition, in the model, and was 
temperature controlled at close to saturated conditions, on 


the actual column. 


eel.45 Variable Holdup 


Several correlations for flow over weirs and bubble 
cap tray holdups have been presented in the literature as 
discussed in Chapter 2. As will be shown in Chapter 5, com- 
parison of predicted holdups with experimentally measured 
holdups from the distillation column used in this study show 
Significant differences, with the measured values being lower 


in magnitude by as much as one half. 


Attempts were made to correlate the experimentally de- 
termined holdups at different steady state conditions with 
liquid and vapor flow rates predicted by the constant holdup 
mathematical model developed in this work, but due to the 
Beatrer Of gata s( seeeAppendix (fF). Figure (5.10) the corre— 
lations were not of value. However, the data did show that 
when plotting measured tray holdups against column feed rate, 
a large fraction of the total column holdup occurred below 
the feed tray, and furthermore, when the feed rate was in- 


creased, this proportion increased. 


For these reasons, correlations presented in the litera- 


ture were modified to give holdups in the same range as 
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measured values. The equations were of the form 
Wwe Rae tia OV 2 et UC 
n nen nian n 


and the modifications were to the magnitude of the parameters 


Bee a and. C>.. 
n n n 


Three correlations were used in simulations in this 
work: in two of the three, the parameters ans oon and Cc, were 
constant from tray to tray and were not changed with time. 
The third correlation (Francis Weir) used parameters that 


varied with density changes. 


Two alternate methods of allowing for variable tray 
holdup were attempted in this work, and are described in the 


Peroblem Formulation” section, which follows. 


3.2 PROBLEM FORMULATION 


Two solution schemes to solve the 4n simultaneous equa- 
tions were formulated. As previously mentioned, 4n unknowns, 
n compositions, n liquid flow rates, n vapor flow rates, and 
m itray holdups *needvto be “specified to completely define the 
column dynamic state at any time. These variables must be 
determined in some way from the composition balance equations, 
total mass balance equations, energy balance equations and 
holdup correlations. The sequence in which these equations 


are solved together to define the column's state, is termed 
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the “problem formulation" in this work. 


Se2e. ) substitution Method 


Seaet.t Description 


In the substitution method, only the composition balances 
were integrated; the other differential equations were 
treated as algebraic equations with their differential terms 
calculated as finite differences. Composition balance inte- 
gration yielded individual tray compositions. The coefficients 
of the liquid and vapor composition terms were evaluated at 


the beginning of the time step. between integration itera- 


tions. The equations integrated were arranged as follows. 
dx 
sien ( x Seley ct pV 4 re es 
at WT, i) eet Meer n=l" n=-1 n Fn 
qwT. 
= = ee Eyes 
Vn¥n a dt one ( ) 


To solve for the liquid and vapor flow rates, the en- 


thalpy balances and total balances were solved simultaneously. 


Total mass balances are given by equation (3.1) 
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where the derivative term is evaluated from holdups of the 


previous time step, from a finite difference representation. 
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Sek eh KL 
= (WT - WIS”) /DT (3.19) 


The energy balances were rearranged to the following 


form 
dh awt 
ee eae anne 
ee a eras - Hehe a es - eae eee (3.20) 


with the enthalpy difference calculated from the finite dif- 


ference representation. 


K 
dh 

ge) ag Re eae 
Spee, ay esas ty Dan AeM (3.21) 
The set of equations (3.1), (3.20) were solved by sub- 


pebcut son fOr the, liquid: and! vapor flow rates. Substituting 
Souatrons (3.1) into-equations (3.20), starting at stage. 9, 
eliminates all the liquid flows, except the known reflux 

flow (Lig). so the resulting equations, in vapor flow only, 
may again be substituted into each other. Upon this, one 
equation in one vapor flow is obtained, and this equation is 
solved. Back substitution then gives all liquid and vapor 
flows for all stages. A detailed step by step procedure with 


detailed equations is included in Appendix (C). 


Accompanying the change in liquid and vapor flows is a 
change in tray holdup, established from equations (3.4). As 


liquid accumulating on the tray has the enthalpy of the 
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liquid leaving the tray, all accumulation comes from the 
existing liquid for the enthalpy balance to be satisfied. 
This correction must be made using an iterative technique; 
holdups are calculated based on new liquid and vapor flows, 
accumulation rates are calculated with equation (3.21), and 
Premexi ting ibigqguid wilow isscorrected for accumulation... A 
holdup is again calculated based on the corrected flow, and 
the process is continued until the difference between succes- 


Sive holdup calculations falls within some error criterion. 


Underlying the entire solution method is the condition 
that the holdup derivative, upon which the calculations of 
current liquid and vapor flows and holdups are based, is cal- 
culated from the previous time step (equation 3.21). Since 
the holdup derivative is defined based on holdup values at 
the K and K-l time instants, and is used to calculate flow 
rates and holdups at time step K+l, this naturally leads to 
some inaccuracies. This is overcome by using the matrix in- 
version method, which will be considered later in this 


chapter. 


3.2.1.2 Stability 


Numerical stability of this solution scheme can be 
divided into two categories, namely: the stability of the 
integration and the stability of the coefficient determina- 
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Once the problem has been formulated into a format such 


as 
reseed G.920) 


there are many standard methods to evaluate stability and 
obtain a stable integration for a variety of A Matrices =the 
primary obstacle in modelling systems is to formulate the 
problem into a standard format. For distillation systems, 
this matrix is near ideal, strongly diagonally dominant, with 
well conditioned eigenvalues, such that a simple integration 
scheme such as Euler or Modified Euler will suffice fora 
stable solution. A detailed description of the components of 
Eis Matrix, which iS important in the matrix inversion method 
SersoOlution, can be found in section 3.2.2.2. “Furthermore, 
this matrix represents the composition system of the distil- 
lation column, whose eigenvalues are inherently smaller than 
those of the hydraulic system, which is essentially the co- 


efficient determination system. 


Although the solution for the A matrix was obtained by 
solving a series of algebraic equations, in essence, inte- 
grations were performed here also. However, the equations 


could not be formulated into a standard form expressed as 


K+1 K+1 K Oey 


so only an ad hoc stability analysis could be done. The 


Sie) 


dependence of the new value of A upon the old value comes 
through the enthalpy, holdup and composition derivative terms 
in the equations. If it were not for the inclusion of the 
holdup derivative, which is independent from the concentra- 
tion derivative, stability of the coefficient determination 
step would be identical to that of the integration step, as 
all differential terms would be linear functions of the con- 


centration derivative. 


The solution with variable holdup resulted in a reduc- 
tion in stability from the constant holdup case. It was 
found, that for the variable holdup case, a time step approxi- 
mately 60 times smaller than for the constant holdup case was 


necessary for a stable solution. 


The change in time step size was considered due to two 
effects: firstly, as previously mentioned, a relatively poor 
estimation of the holdup derivative was used, as only old 
holdup values were available for its calculation. The error 
incurred was necessarily reduced by using a smaller step 
size. Secondly, the hydraulic system, to which the tray 
holdups are directly related is a much faster system, as com- 
pared to the concentration system, with relatively shorter 
time constants. In light of this, smaller integration time 


steps would be expected. 


The necessity of employing a very small integration 
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interval provided the incentive for restructuring the formu- 
lation leading to the matrix inversion method that is con- 


Sidered next. 


3.2.2 The Matrix Inversion Method 
9.2.2.1 Description 


In order to reduce solution time due to the small time 
Mncrement,! it: was decided that a better estimate of the ave- 
rage tray holdup rate of change, over the integration inter- 
val, was required. An exact expression was available in the 


Eotal+ mass balance (equation 3.1). 
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This set of equations was integrated directly in an integra- 
Bron step, along! with the component balances (3.2), to pro- 
vide trajectories for tray compositions and tray holdups. 

The vapor and liguid flow rates were again calculated ina 
coefficient determination step. The energy balance equations 
and holdup correlations were rearranged and used for this 


purpose. The energy balance (3.3) may be rewritten 
é 


qwrt dh 
ee - - V H 
a dt 2 on ae sere rul ees nin 
= h - (ctu. 2) 
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Replacing the rate of change of holdup using equation (3.1) 


and rearranging yields the equations (3.25) in the required 
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As opposed to the substitution ee this method uses 
a current value of the holdup derivative based on current 
Pomues OL liquid and vapor flows. By substitution into both 
the enthalpy balance and composition balance, the holdup de- 
rivative and liquid flow terms for a tray were eliminated 
from the corresponding equation (see equation 3.25). A de- 


tailed derivation is included in Appendix (C). 


The holdup correlations were also rearranged (equation 
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When equations (3.25) are combined with equations (3.4), a 


matrix equation is the result, 


(Sane oy 


where FT is the vector comprised of the left hand sides of 
the above mentioned equations, and TT is a matrix of co- 
efficients which are enthalpy differences and holdup corre- 


lation parameters. 
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TO solve: for Liquid and vapor flow rates, the vector FT 
was updated for changes in holdup, and for changes in en- 
Eialpy accumulation. Inversion of matrix TT, and taking the 
product of the inverse and the FT vector permitted calcula- 


mon Of the mass flow rates from 


PT (3h25) 
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Thus, in one time step, new liquid and vapor flows were 
calculated, based on current compositions and holdups. Then 
via integration, new compositions and holdups were determined 


from the new mass flow rates, that is 
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The integration step included the composition integra- 


tion, according to the following expression derived in 


Appendix (C). 
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To calculate new holdups, the first order Euler inte- 


gration was used 
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3.2.2.2 Stability 


Due to the nature-of the equations, rigorous stability 
analyses could not be made. For the integration step, but 
not for the coefficient determination, the equation's format 


was 
on = eh (2 tes) 


to which standard methods could be applied. For the inte- 
gration step, the A matrix was formed as a tridiagonal matrix, 
its derivation being presented in Appendix (C). However, as 
will be shown later, it was found that the stability of the 
integration step, as in the substitution method, was not the 
controlling factor. Rather, the coefficient determination 
step, viewed as representing the hydraulics, limited the size 


of the integration interval. 


It was found that using the identical holdup correlation 
used in the substitution method, that for the Matrix Inver- 
sion method, a time step of approximately 3 to 4 times larger 
could be used to obtain a stable solution. However, this 
time step was very dependent upon the liquid coefficient "a," 
mm the holdup corretation (3.4), with a larger "a" permitting 


a larger time step to be used. The reason for this was evi- 


dent upon inspection of the equations. 


Wire ee far tb V 
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A small change in holdup would change the liquid or 
Vapor flows in order to satisfy the holdup correlation. As 
the enthalpy balance is also satisfied, and since liquid en- 
thalpies are lower than those of vapor by approximately a fac- 
tor of 10, the liquid flow rate would absorb most of the change. 


So, Or example, jiré "a were small, a change on the left hand 
Side of equation (3.4) would cause a relatively large change 
in liquid flow. This change, over an integration interval, 


must for stability be limited by using a smaller time step. 


All of the results and comparisons shown in Chapter 5 


were obtained by the more stable matrix inversion method. 


3.3 SIMPLIFYING ASSUMPTIONS 


To compare the response of the rigorous model with dyna- 
mic responses of models based on additional assumptions, 
changes to the detailed model as described in the following 
sections, were introduced. The matrix inversion method of 
solution, suitably modified for the less rigorous models, was 
employed, as it is simpler, more straightforward and faster 


than the substitution method. 
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and the total mass balances, equation (3.1), formerly diffe- 


rential equations become algebraic 


7 SRU Mira VEER. Be Suck Poe eek (ey32) 


n+l 


Two cases of constant holdup may be assumed, constant 
volume holdup and constant mass holdup. For each, the right 
evae OT ‘equation (3.32) assumes a different value. In the 
former case of constant volume holdup, mass accumulation will 
be due to variation in composition and thus a density change. 
As density is a linear function of composition, the following 
expression holds. 

dx 

K =c == (3313.3) 

For the assumption of constant mass holdup, it follows that 


K = 0 (S334) 


In this work, only the case of constant mass holdup was 
considered. Changes to implement this assumption were made 
to both the integration and coefficient determination sys- 
tems. The integration of the total mass balance equations 
in the integration section was dropped. These equations were 
used in place of the holdup correlation equations in the 
solution of liquid and vapor flows. The matrix "TT" and vec- 
tor "FT" from equation (3.26) were changed accordingly. 

Apart from these changes, the solution in each time step was 


Carried out as for the variable holdup model. 


ve PL i Ad 
18 DLE ReMae os 
; : - - = 


* f " 


40 


3.3.2 Constant Liquid and Vapor Rates 


The variable holdup model may be further simplified by 
assuming that the liquid and vapor rates do not deviate 
greatly from their steady state values for feed, reflux or 
steam flow disturbances of moderate magnitude. As holdup is 
a tunction!of’liquid and vapor rates, constant holdup’ must 


also be assumed. 


These simplifications, which lead to a reduction of 
order from forty states to ten, reduces computation time, as 
only the integration of equations (3.31) remains. Also, 
matrix A remains constant except for changes in the distur- 


bance variables. 


wo. oe Constant Molar Overflow 


If the molar heats of vaporization of the two components 
are nearly equal, a mole of vapor condensing will cause a 
mole of liquid to vaporize. If heat losses from the column 
are neglected, then flow rates, both vapor and liquid, will 
be uniform throughout the stripping section and again through- 
out the rectifying section. The difference between the flows 
in the rectifying section and stripping section will be 
solely due to the feed flow, which is introduced between the 
two sections. This assumption eliminates the need for an 


energy balance for each tray. 


4 
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Vapor dynamics were taken to be instantaneous resulting 
from the assumption of constant vapor holdup. Liquid dyna- 
mics, however, were handled by including the holdup corre- 


factions in-the solution scheme. 


Modifications to the matrix inversion solution method 
were as follows. In the energy balances, the vapor and liquid 
enthalpies, the coefficients for the vapor and liquid flow 
terms, respectively, were taken to be fixed values, for a 
composition of 50 percent MeOH. This resulted in using con- 
Stant molar heats of vaporization from tray to tray. Heat 
input to the reboiler was calculated as before, and largely 
determines the vapor rate throughout the column. The holdup 
correlations, total and component balances were included in 


the solution scheme as before. 


B.o.4 Adiabatic Operation 


in dist. Viation columns, sthere.is heataloss -from athe 
column to the atmosphere by radiation and convection. A 
further simplifying assumption that may be employed, is the 
neglect of heat loss from the column. Mathematically, this 
results in dropping the heat loss (Q term) from the energy 


balance equation (3.3) yielding 
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dh qwrt 
en dt pga tet at i ona ve Ln 
= Ven - Veqeaci ~ Sha + F pen (34.3 5)) 


For columns housed within buildings such as the one at 
the University of Alberta, heat loss is largely a function 
of column wall temperature. A significant portion of the 
total heat loss from the column was attributed to the gas 
chromatograph sampling system, where a long length of metal 
tubing, through which hot reboiler liquid was cycled, was 
exposed. This was a Loss directly from the reboiler, amounting 
to sae three to four percent of the total heat in- 
put to the column. Heat loss from the remainder of the 
column made up the balance, which amounted to about three to 


four percent of the heat input. 


In the simulations, an identical heat loss for each tray 
was employed, with the loss selected so that the terminal 
compositions matched the experimental values at the initial 
steady state. When adiabatic operation was assumed, initial 
steady state conditions could not be matched, indicating the 


severity of this assumption. 


3.3.5 Variable Heat Loss and Tray Efficiency 


As a distillation column experiences upsets affecting 


liquid and vapor compositions on the stages, the temperatures 
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of the stages and thus individual tray heat losses are also 
affected. In the previous section, the assumption of adia- 
batic operation was described. An alternate assumption, 
rather than consider the heat loss to be constant, is to 


@etOw=f0ra heat. loss =that varies with operating conditions. 


The constant heat loss assumption would mathematically 
imply that the Q term in the energy balance, equation (3.3) 
is taken at a constant value over the course of the column 


response to a disturbance, that is 


Rather than use a constant heat loss value, heat loss 
was varied directly with composition. Equation (3.37) was 
used to vary the heat loss with composition changes according 


to 


: PS 
are ae Glaeis 
x 


where xy is a composition of a chosen stage in the column at 
time t.- In this work, the» chosen stages were taken as the 
feed tray, the reboiler and the condenser. The heat loss at 


the initial and final steady state conditions were chosen to 


match the experimental steady state conditions. 


Tray efficiency, like heat loss, was also allowed to 


vary with composition changes. Values were selected to match 
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experimental and simulated initial and final steady state 
terminal compositions. Efficiency was varied, in a manner 


directly amalagous to that for heat loss, according to 
X pig 
Bye 
ai oar a vee 


This equation was used for both the stripping and rectifying 
sections, as efficiencies were taken to be uniform through 


a section, but different for the two sections. 


3.4 NUMERICAL METHODS 


Several different methods were used to integrate the 
composition balance equations. These equations were inte- 
grated independently in the substitution method, and as a 


matrix system in the matrix inversion method. 


4 4) The Substitution, Method 


Equation (3.18) was expressed as below for the integra- 
tion purposes for any plate other than the feed plate and 
with no sidestream withdrawal. 
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For the feed plate (stage 5) the equation became 


dx. l 
at Wi, eee tests GYs VY ava 
awT. 
oP PX, = ac ar C3 4.0:) 


The equation for the reboiler (stage 1) became 


X57 Lj Xp ~ Vay ae (3.41) 


The equation for the condenser (stage 10) is 


= X19)? (3.42) 


At any time instant, K, the general vector-matrix re- 


presentation may be expressed as 


K 
ax 
—n K 
a5 S343 
dt tn ( ) 


The Euler method is summarized by equations (3.44) and 


45) 
K 
“En Car = a, 0h (3.44) 
rearranging 
gh Tess ei pa ER) (3.45) 
In this simple integration scheme, equation (3.44), an 


elemental but poor approximation of the actual derivative, 


is employed. In calculating a transient, an error is incurred 
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at each time step, each error being propagated over the rest 
of the calculation. The modified Euler integration method 
offers a slight improvement over the simple Euler integration 
in the estimation of the derivative, and allows larger inte- 


gration time intervals to be used, at the expense of more 


computation. 
axk*? ax* axktt 
fia ge deen eat sop 


The derivatives on the right hand side of equation 
(3.46) evaluated at times K and K+l are calculated using the 
Euler technique after which an average derivative over the 
time interval is calculated and used for the integration 


interval. 


As a further attempt to better estimate the value of the 
composition derivative, third and fourth order Runge-Kutta 


techniques were used; they are outlined below. Referring to 


the £ function described in equation (3.43), let 
a S (3.47) 
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then 


il 
— + — 
*n sel 6 ' din i aI Fee ui ahs) 
Note that: equation (3.47) does not include time as an 
independent variable, and thus time does not enter into the 
fourth order Runge-Kutta solution described above, as a 


Variable. 


The more complex integration schemes such as the Runge- 
Kutta methods allowed the integration time interval to be 
increased, but involved increased computation. But, as men- 
pronediin sections 3 32. 1..:2..and 3.2.2.2, the integration steps 
in the model solutions were not integration interval size 
Mimlcang; rather, at was the step of calculating vapor and 
liquid flow rates and liquid holdups that forced the use of 
the small integration step size. The integration method, 
therefore was not crucial, and it was found that the simple 


Euler method served as well as the more complex techniques. 


3.4.2 Matrix Inversion Method 


As outlined in section 3.2.2, the matrix inversion 
formulation dealt with the integration as a system of equa- 


tions expressed in matrix form as 
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Use of the Euler method for integrating the system equa- 
tions expressed in matrix form as given above yields 
Gh = x) /pt = AXx* 


ee DT A) x (3154 8) 


Also employed was the modified Euler method leading to 


ax aie + aS 

ae ae es) 

at 8 2 

eel tena nett Ax (3.49) 


The derivation is given in Appendix C. 
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CHAPTER 4 


SIMULATION RESULTS AND DISCUSSION 


4.1 RESULTS - INTRODUCTION 


Results of distillation column simulations demonstrating 
the effects of various simplifications and assumptions are 
presented. The basic model for the ten stage binary column 
was chosen to the 30th order, in which the total mass, energy 
and component balances were solved simultaneously. To allow 
Boreatvariable holdup, introduction of holdup correlation 
equations into the solution aay cad in a 40th order model. 
Results for three cases of this model are presented in this 
chapter. The responses of various different models of lower 
complexity relative to the basic model are also presented for 
comparison. These are, in order of increasing complexity, 
the constant flow rate model, the constant molal overflow 
model, and the constant efficiency and heat loss model. In 
addition to the simulations for the basic model, for which 
the parameters of tray efficiency and heat loss are taken as 
functions of distillate composition, results are also pre- 
sented for the case where the parameters are functions of feed 
tray composition and bottoms product composition. Finally, 
simulated results obtained including as heat sinks the tray 


metal and glass in the basic model are also presented. 


The figures displayed in this chapter enable comparisons 
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of the effects of various assumptions upon distillation 
column dynamic responses. These responses further demonstrate 
the robust nature of the solution scheme outlined in Chapter 


3, and its applicability under a variety of conditions. 


A complete summary of the simulation results is presented 
in Table 4.1, along with an explanation of the codes used in 


Peqgures 4,1,.to 4.12 inclusive. 


Results are presented for step disturbances of -20 percent 
in feed and reflux flow rate and +10 percent in steam flow rate. 
Only the disturbances of reflux and steam were chosen with which 
to do complete comparisons for reasons of time and expense. 
These disturbances were chosen since it was found, as will be 
shown in Chapter 5, that these two disturbances exhibited the 
largest discrepancies between experimental and simulated res- 
ponses. Furthermore, in many previous modelling studies only dis- 


turbances in feed composition or feed flow were considered. 


The model using constant holdup, variable efficiency and 
heat loss (the 04 series, Table 4.1) is considered as the 
reference (base) model for purposes of comparison to predicted 
responses. This is a more general model than that of Huckaba 
et ae which is considered a general model, representing data 
well. Responses calculated from models of both increasing or 
decreasing complexity are compared with the base case res- 


ponses. 
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TABLE 4.1 


SUMMARY OF SIMULATIONS 


I Reflux Disturbance 
(Minus 20 percent) 


Run Code Assumptions Figure 
RO1M Ly sie 4.4 
RO2M Pipe paile ay it | 4.4 
RO3M Soe case 4.4 
RO4M oN Meee Be 4.4 
RO5M eo ilD 4.6 
RO6M Syeb2ey Ls 4.6 
RO7M 4,7 4.9 
RO8M 4 LZ 4.6 
ROOM 4,9 4.10 
R1LOM Ay.S 4.10 
R11M A lduytlS 4.9 


III Reflux Disturbance 
(Plus 20 percent) 


Run Code Assumptions Figure 
RO3P Sess 4.2 
RO4P Aides a2 

V Steam Disturbance 

(Minus 10 percent) 

Run Code Assumptions Figure 
SO3M Pt es 1) 4.3 
SO4M Re eas LS) Aes 


II Steam Disturbance 
(Plus 10 percent) 


Run Code Assumptions 
SO1P Me lee S 
SO2P PAS, COS RIS) 
SO3P Bie 2ids 
SO4P AL 2S 
SOSP ogee. 

SO6P Gipnire 
SO7P Are 
SO8P 4,12 
SO9P 4,9 
S1OP 4,8 
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IV Feed Disturbance (Flow) 


(Plus 20 percent) 


Run Code Assumptions 
FO3P EA ae Ee. 
FO4P AD Mie pues 


Figure 
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TABLE 4.2 


LIST OF ASSUMPTIONS 


constant molal overflow 
constant liquid and vapor rates 
constant efficiency and heat loss 


variable efficiency and heat loss: based on top 
composition 


variable efficiency and heat loss: based on bottom 
composition 


variable efficiency and heat loss: based on feed 
tray composition 


variable holdup: arbitrary correlation 
variable holdup: RADEMAKER correlation 
variable holdup: FRANCIS WIER correlation 
subcooled reflux 

subcooled feed 

constant tray liquid mass holdup 


negligible heat capacity of tray glass, metal 
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4.2 SIMPLIFYING ASSUMPTIONS 


week constant Breiclency, Constant Heat Loss 


BseMmentioned in Section 3.4.4, the»parameters* of tray 
efficiency and tray heat loss were used to fit simulated res- 
Ponse@s OO experimental data. Figures 4.1 to 4.3 present a 
comparison of predicted responses obtained by assuming con- 
Stant tray heat loss and constant tray efficiency (the 03 
series) with responses obtained by allowing for changes in 
heat loss and efficiency (the 04 series) for feed, steam and 
reflux disturbances. It can be seen that in the constant 
parameter case, the responses exhibit the same general dynamic 
behaviour as in the variable parameter case, but a compari- 
son of response gains shows far greater gain values for the 


case of constant parameters. 


This result is to be expected. Heat loss, which varies 
with tray temperature (or composition), varies in such a 
manner to reduce the effects of a disturbance; as a distur- 
bance causes column temperatures to rise, heat loss rises to 
compensate, tending to reduce the effects of the disturbance. 
The same is true for disturbances causing temperature reduc- 
tions. However, correlations relating efficiency to column 
parameters indicate a dependence upon flows and holdup; 
changes in these variables possibly do not have a damping 
effect. Heat loss, however, has been found to have much more 


effect on column operation than efficiency, and it is mainly 
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FIGURE 4.1 +20% Feed: Effects of the assumption of constant 


heat loss and efficiency. 
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FIGURE 4.2 +20% Reflux: Effects of the assumption of constant 
heat loss and efficiency 
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the damping effect of the varying heat loss that accounts 


for the gain differences observed. 


That the dynamics are unaffected may also be expected. 
Heat loss and efficiency have negligible dynamics of their 
own, but vary directly with column parameters such as flow 
and temperature. In all cases, the temperature (or composi- 
tion) dynamics are slower, and thus largely determine res- 


monse rates. 


aie Constant Flow Profile 


The assumption of constant flow profile is a severe 
assumption, causing liquid and vapor rates to be, not time 
dependent, but different from tray to tray. The effect of 
this assumption may be observed for the product composition 
responses for reflux and steam disturbances shown in Figures 


4.4 and 4.5 respectively, as the 02 series. 


The constant flow profile model predicts responses with 
lower gain and slower dynamics for both steam and reflux up- 
sets. Such disturbances to a column cause compounded effects 
Of flow, composition and energy. Since in the solution 
scheme of this model, the tray energy balance and total tray 
mass balance equations do not need to be solved, the resulting 
effect is to reduce the magnitude of the simulated column 


upset as shown. The error due to this simplification appears 
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FIGURE 4.4 -20% Reflux: Effects of the simplifying assumptions 
of constant mass flow rate, equal molal overflow 
and constant heat loss and efficiency. 
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FIGURE 4.5 +10% Steam: Effects of the simplifying assumptions 


of constant mass flow rate, equal molal overflow 
and constant heat loss and efficiency. 
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to be much more pronounced in the bottoms composition res- 
ponse than for the distillate response. Although a very 
restrictive assumption, the error incurred by the use of 

this assumption may be acceptable for use in distillate con- 
trol studies, where compositions do not deviate far from the 
initial steady state. The savings in computing time may out- 
weigh the lack of accuracy in representing column nonlineari- 


ties. 


ees 6hCONnStant Molal Overflow 


The constant molal overflow assumption is a severe as- 
sumption if the heats of vaporization are much different. 
For the methanol-water system, the heat of vaporization of 
water is 1.8 times the heat of vaporization of methanol. By 
assuming an average value of the heat of vaporization, simu- 
lations predict more vapor flow in the rectifying section of 
the column and less vapor flow in the stripping section, than 


would actually occur in the pilot plant column. 


Phe errors due to this simplification may be seen In: the 
01 series in Figures 4.4 and 4.5 which show the predicted 
top composition responses for both reflux and steam distur- 
bances. In the case of top composition, deviations from 
steady state are greater than for the case of constant liquid 
and vapor flow, and are close to the cases of constant 


efficiency and heat loss. In the constant molal overflow 
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case, the flow profile throughout the column will change 
directly with composition changes only; enthalpy effects, 
occurring along with composition effects in an actual column, 


are averaged out. 


The bottoms composition responses in the case of con- 
stant molal overflow are, as before, close to the other si- 
mulated responses. This is occurring for sensitivity rea- 


sons, mentioned in the following section. 


4.2.4 Heat Loss, Efficiency Functions 


As previously mentioned in section 3.3.5, the parameters 
of tray heat loss and tray efficiency were varied with the 
composition of a single stage in the column, during the 
duration of the response. In actual operation, these para- 
meters likely vary individually with tray compositions and 
Elow rates. For simplification, however, these parameters 
of heat loss and efficiency were taken as functions of a com- 
position at one point in the column. This point was changed 
to three different locations in three separate simulations: 
the condenser, the feed tray, and the reboiler. The effects 


of these changes may be observed in Figures 4.6 and 4.7. 


Fitting inaccuracies are shown in the lack of agreement 
of compositions at the final steady state. Moving the refe- 


rence location for the heat loss and efficiency functions 
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heat loss. 
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will only change response dynamics and not gains. For ini- 
tial and final steady states of simulations of different 
cases of the same model to agree, identical initial and 
final efficiency and heat loss values must be used for all 
cases. The parameters, P, of heat loss and efficiency, were 
varied as expressed by Equation 4.1. 

(xX, - x,) 


P, = P. + ——— (P. - P.) (ao 1) 
ce 1 (X- - x.) 1 


Although the desired final efficiencies and heat loss 
were chosen to match the final experimental steady state 
compositions, and were constant for all three simulations, 
calculated parameters at the final steady state (Pi, t = t.) 
did not exactly match the desired final parameters, and 


caused the discrepancies. 


A trial and error method was used to obtain final steady 
state efficiency and heat loss terms that would enable 
matching the simulated terminal compositions to experimental 
ones. A summary of efficiencies and heat loss values that 
permitted matching for each disturbance is presented in Table 


ees 


However, it is still valid to make qualitative observa- 
tions as to how dynamics are altered by changing the reference 
location. In all cases, the dynamic behavior of the bottom 


composition responses are almost identical to the base case 
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TABLE 4.3 


SIMULATION HEAT LOSS AND EFFICIENCY 


INITIAL HEAT LOSSES INITIAL EFFICIENCIES 

Trays 400 joules/sec Sori pping 5S 

Reboiler 1000 joules/sec Recerlrying ssw 

DISTURBANCE FINAL HEAT LOSS (j/s) FINAL EFFICIENCIES 
Trays Reboiler Stripping Rectifying 

+10% steam 700 1000 95 .90 

-10% steam 200 800 mesiis ee 

+20% reflux 190 1000 .80 TO 

-20% reflux 600 1000 viet ahe 

+20% feed 250 800 93 Swe: 

-20% feed 650 1000 293 .o5 

(04) series). Disturbances examined forced the bottoms com- 


position into an insensitive region, at approximately 0.5% 
methanol, in the pinch region of the equilibrium curve. As a 
result, the gains were minimally affected. In addition, the 
relatively high circulation rate in the stripping section and 
the large holdup in the reboiler combine to make the reboiler 
composition dynamics dominant. These dynamics largely depend 
upon reboiler holdup so not much change from model to model 


would be expected. 


The effect of a changing reference composition location 
does however, significantly affect the top composition dynamics. 


The three different reference locations for varying heat loss 
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and efficiency exhibit three distinct effects that are common 
to both steam and reflux disturbances. With heat loss and 
efficiency varying with feed tray composition (06 series), dis- 
tillate composition responds most rapidly, since the feed tray 
reacts quickly to disturbances, the feed tray composition being 
in a very sensitive region of the equilibrium curve. With heat 
loss and efficiency varying with bottoms composition (05 series) 
distillate composition dynamics show the dynamics of two ef- 
fects; the flow of energy dynamics cause the rapid reaction of 
the bottom product composition to the flow disturbance. As the 
bottom product composition more closely reaches its final steady 
state, the composition effects working their way slowly through 
the column, react more slowly at the condenser, eventually 


trimming all compositions to their final steady state values. 


The smoothest of the three responses, as is to be expected, 
Securs for the case of heat loss and efficiency varying with 
distillate composition (04 series) which seems to have the 


slowest response to disturbance in all cases. 


The top product composition responses for negative reflux 
flow disturbances, however, exhibit an unexpected nonlinear 
type of behavior that is not typical of the experimental res- 
ponses (as seen in Chapter 5). It is considered that by letting 
heat loss and efficiency vary with a certain composition in the 
distillation column, the dynamics of that particular stage are 


directly introduced and amplified in the dynamics of other 
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compositions in the column, as the simulated responses show. 
In addition, the assumption of uniform heat loss and uniform 
efficiency in the rectifying and stripping sections cause the 
Simulated column to adopt unnatural flow and composition pro- 


files, perhaps also contributing to the different curve shapes. 


A model of further complexity could be developed, with 
heat loss and efficiency varying individually with tray com- 
positions or temperatures. The fitting involved in such a 
case, however, would be a very complex and time consuming 


procedure. 


fee. oo thermal Capacitysof Tray Metal sand Glass 


The possibility that heat effects due to the energy 
storage in the tray metal and glass could be significant, was 
examined. The modifications to the model to include this 
effect involved the addition of one extra term in the energy 


balance equation, the last term shown in Equation (4.2). 
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The plots of predicted composition responses allowing 
fOr this effect are_shown in Figures 4.6 and 4.8 (the 08 
series). Comparison of the responses in these figures with 
those of the base case show such effects to be negligible, 
for the pilot plant column used in this study. Changes in 
the responses may be observed as slight smoothing, evident 
upon close examination. Also evident are the marginally 
peower Overall dynamics.. Both of these effects are expected, 
an additional heat reservoir having a damping effect, and 
may be observed in both distillate and bottoms composition 


responses. 


4.2.6 Variable Holdup 


Responses of models including the assumptions of variable 
holdup are presented in comparison with the base case (the 
04 series) in Figures 4.9 and 4.10 for reflux disturbances 
Figure 4.11 for the steam disturbance. Three different hold- 
up correlations were used. One of those chosen was used by 
Rademaker et Aa cts their work in distillation column model- 
ling; the Francis wier formula for flow over rectangular 


wiers was used in other simulations. The third correlation 
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+10% Steam: Effect of considering tray and glass 
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FIGURE 4.9 -20% Reflux: Effects of the assumption of variable 


holdup according to arbitrary holdup correlation. 
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-20% Reflux: Effects of the assumption of variable 
holdup according to the FRANCIS wier, RADEMAKER 
holdup correlations. 
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FIGURE 4.11 +10% Steam: Effects of the assumption of variable 
holdup according to the arbitrary, FRANCIS wier, 


RADEMAKER holdup correlations. 


ts 


was, as the other two, of the general form presented in 
Equation (4.3), but with arbitrarily chosen parameters ane 


oe rand Cc . 
n n 


wre = at. - SLA oa oo (4.3) 


For the Rademaker correlation, but not for the third 
(arbitrary parameters) holdup correlation, the coefficient 
of the vapor flow, Doe had a larger magnitude than the liquid 
coefficient, an: This was to have an important effect upon 
model trajectories; the basis for the selection of coeffi- 
cient values is presented in section 3.1.4 and in Appendix Cc. 
The vapor coefficient, Doe may be shown to be null in the 
Francis wier correlation, flow over a wier assumed to have 
no dependence upon vapor flow. Table 4.4 presents a summary 


Orc the coefficients of Equation (4.3) used in this study. 


A comparison of the simulations with holdups based on 
mae third correlation (arbitrary parameters, the 07 series) 
with the simulations of the base case shows very little 
effect of the variable holdup assumption upon column dynamics. 
From the bottoms composition response due to a reflux dis- 
turbance, dynamics appear to have been slowed down slightly, 
as is also the case, but more apparent in the steam upset 
simulation (07 vs 04). However, in Figure 4.10, including 
the Rademaker correlation with smaller coefficients in the 


holdup equation resulted in the prediction of responses with 
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TABLE 4.4 


HOLDUP CORRELATION PARAMETERS 


Holdup Series Liquid Vapor Constant 
Correlation Coefficient Coefficient 
a b Cc 
n n n 
- RADEMAKER ET AL 10 Zee -1.54 1500 
FRANCIS? 9 eon Ge 0 vegas 
WIER 
ARBITRARY 7 10 -20 L500 


1 correlation depends upon the following parameters 
tes oud. flow .(g/s) 


0 - tray liquid density (g/1) 


Derivation of the Francis wier expression is presented in 


Appendix F. 


i i 


faster dynamics. In Figure 4.11, the same observation may 
be made with regard to the steam disturbance. Use of the 
Francis Wier formula (the 09 series) resulted in slower dyna- 


mics for both steam and reflux upsets. 


These observations are explained by considering the 
effects of the coefficients iets and ok in Equation). (423) 
upon the holdup and flow rates when a disturbance takes 
place. Flow rates and rates of change of holdup at any time 
are governed by the total mass and enthalpy tray balances, 
while absolute holdups and flow rates are related by Equation 
(4.3). As the enthalpy of tray holdup and liquid leaving 
the tray are identical, a disturbance which forces the hold- 
mp ytoVchange wil) quickly ‘cause the liguid flow to change to 
Satisfy the enthalpy balance. To what extent, and at what 


rate depends upon the holdup correlation coefficients. 


The three sets of equations mentioned above, used in the 
solution of flow rates and holdups, are solved simultaneously 
and it is difficult to describe the propagation mechanism. 
The initial disturbance will cause a change in liquid flow 
from a particular tray, a change in vapor flow, and a change 
in the rate of accumulation of holdup, with the magnitude of 
the changes determined by the values of holdup coefficients 
and the enthalpies throughout the column. To determine the 


exact nature of the individual influences of these parameters 
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would involve a more complete set of simulated LUNs; -LOr 6a 
spectrum of disturbance types and coefficient values. For 
this reason, no quantitative relationships of dynamics and 
coefficient values will be presented here. Rather, some 
examples will show why coefficient sign, magnitude and vapor 


to liquid coefficient—ratio are important. 


Figure 4.11 shows the response of the distillation 
column to a steam disturbance, an upset which initially 
causes mainly an increase in vapor rate. Generally, flow 
changes occur more quickly than holdup changes; because of 
this and the opposite signs of the coefficients of the Rade- 
maker correlation, including the Rademaker correlation in 
the solution speeds up simulated dynamics. In the short 
term, because of slowly reacting holdups, the nature of 
Equation (4.3) would cause a vapor rate increase to predict 
a liquid rate increase, as a result of including this equa- 
tion in the model solution. This results in an increased 
mass circulation rate, spreading the disturbance effects more 


quickly, and speeding up column dynamics. 


Figure 4.10 shows the effects of holdup coefficient 
magnitude upon column response. Large holdup correlation 
coefficients such as used in Francis weir correlation, or 
such as the arbitrarily chosen coefficients ensure holdups 


sensitive to liquid and vapor rates. In a disturbance, 
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holdups change faster and further, absorbing or damping the 
liquid flow changes, and slowing the rate of column response 


(see section 3.2.2.2). 


Figure 4.9 demonstrates the significance of the ratio 
Of liguidecoefficients®tolvapor coefficients: | Fora reflux 
disturbance, almost no effect is shown in column response by 
assuming the arbitrary parameters represent holdup behavior, 
where for the steam a distinct effect is shown. In both 
types of disturbances there are increases in vapor and liquid 
flow, but in different proportions. The coefficients happened 
mQebe chosen in sucha way, that the vapor and liquid flow 
rate terms in equation (4.3) were not changed by tray holdup 
emanges to affect column dynamics significantly. For other 
disturbances, however, these coefficients had a pronounced 


effect upon column dynamics. 


The Francis weir correlation predicts a zero dependence 
Of liquid holdup upon vapor flow, and thus vapor effects are 
eliminated. Furthermore, the correlation also predicts abso- 
lute holdup values approximately double the average measured 
values, thus giving tray time constants that are double 
actual values. This results in a severe damping of liquid 
flow and composition dynamics as can be seen in Figures 4.10 


and 4,121. 
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coefficient magnitude, sign, and ratio are important in 
determining simulated column dynamics. The speed of the 
dynamics may be altered by choosing different sets of holdup 


correlation coefficients. 


In conclusion, the effects of allowing a variable volume 
holdup in the model will depend upon the holdup correlation 
chosen. Simulation results in this work have shown these 
effects to cause an increase or a decrease in the distilla- 
tion column models rate of response to a disturbance. The 
results also show that the oan scheme of the nonlinear 
distillation column model developed in this work is stable 


enaevalia for a variety of holdup correlations. 


4.2.7 Enthalpy Reductions 


The feed and reflux enthalpies were reduced by 5 kJ/kg 
each in independent runs, to determine the effect upon column 
dynamics. The results of these changes would be to raise 
the initial steady state methanol compositions on all trays. 
The results, presented in Figure 4.12 for steam disturbances, 


indicate little change in dynamics. 
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CHAPTER 5 


MODEL VERIFICATION - RESULTS AND DISCUSSION 


5.1 EQUIPMENT 


ae) ‘The. -Colwumn 


In order to verify composition transients predicted by 
the mathematical model, open loop tests were conducted on a 
pilot scale distillation column separating a methanol-water 
feed. The column configuration, shown schematically in Figure 
9.1, is essentially the same as that used by recent control 
study workers (28, 44, 47 ) at the University of 
Alberta, except for changes in bottoms composition measure- 


ment system”as described “in section’5.1.2. 


The 8 tray, 9 inch diameter column is fully instrumented, 
all pertinent variables being measured or controlled. The 
walls are glass enabling viewing and estimation of froth 
heights. The stainless steel bubble cap trays were built at 
the University of Alberta, to accommodate thermocouples and 
fluid circulation for composition measurement purposes. The 
column is equipped with a total condenser and thermosyphon 
reboiler. Constant liquid level in the condenser and reboiler 
is maintained by adjusting product flow rates, distillate at 
the top and bottoms from the reboiler. The feed rate, reflux 
rate and steam rate, which are the major manipulated and/or 


disturbance variables are measured and controlled at their 
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set points. Cooling water flow rate to the condenser is em- 
ployed to regulate tower pressure. Also controlled are the 
inlet temperatures of the feed and reflux streams. Plate 
liquid temperatures, cooling water temperatures, steam and 
condensate temperatures, product temperatures and several 
temperatures in the reboiler are monitored. The top product 
composition, high in methanol content was measured by a capa- 
citance probe, and was controlled by an electronic controller 
cascaded to reflux flow rate. An IBM 1800 digital computer 
with a Direct Digital Control package provided control and 
data acquisition, keeping historical records of important 


measurements. 


An online gas chromatograph, used to measure the bottoms 
composition at 4 minute intervals, was driven by the 1800 
computer. It was equipped with an automatic liquid sampling 
valve, a heater, and a special liquid circulation loop. A 


detailed description is described in Appendix E. 


The feed solution was continuously recycled from a large 
storage tank, containing a sufficient quantity to prevent 


feed composition drift during the course of a run. 


5.1.2 Bottoms Composition Measurement 


As methanol-water mixtures high in water content can not 


be analysed reliably by a capacitance system, an alternate 
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method was found. This consisted of a Hewlett Packard model 
5722A chromatograph with an automatic liquid sampling valve | 
and a product circulation system, as shown schematically in 


Pagure 5.2. 


The gas chromatograph is a discrete analysis system; a 
sample is injected into a separation column at periodic in- 
tervals. A carrier gas draws the sample in gaseous form over 
a bed of porous material, having an adsorption selectivity to 
the components of the sample. As a result, the components are 
physically separated as the gas passes through the bed, and 
pass out of the bed as peaks of concentration, normally dis- 
tributed about a characteristic elution time. These impuri- 
ties in the carrier gas cause a change in heat transfer pro- 
perties of the combined gas, and passing the mixture over 
heated metal filaments will cause a temperature rise or fall 
Of the filaments, altering their electrical conductivity. 

The change in filament resistance is converted to a change 

in voltage and from the voltage signal, ratios of peak sizes 
determine relative quantities of components in the sample. 

A sample voltage signal is presented in Figure 5.3 for a sample 
containing approximately 5 percent methanol. The upper plot 

of Figure 5.4 shows only the methanol peak of the same 


sample, showing an absence of noise. 
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For a reliable analysis, it is necessary to have good 
separation of the components in the column. This necessitates 
careful adjustment of gas chromatograph operating conditions, 
column temperature and carrier gas flow rate. The injection 
of sample into the separation column is also very critical as 
the sample should come into contact with the column packing 
as a slug, to prevent tailing of peaks or sliding of one peak 
into another. Also, the sample must be in vapor form at that 
point. Proper sampling requires that a sample of uniform 
Size from the circulating product stream be flashed quickly, 
and introduced as a slug to the separation column. Only one 
of three sampling valves used in the course of this study 
satisfied these requirements. A more thorough account of 


these valves is presented later in this section. 


The sample volume of the first valve used was easily 
enlarged, so that the effect of sample size upon gas 
chromatograph analysis could be determined. The most impor- 
tant and persistent problem encountered was that of repro- 
ducibility and consistency of G.C. measurements. It is known 
that a smaller sample size does give a sharper and more se- 
parated chromatogram, but the voltage signal from the detector 
filaments are lower and, upon reaching the computer analog 
input point, more disturbed by noise. This was balanced by 
the upper limit of column loading; a larger sample gave larger 


peaks but poorer separation. By varying the sample volume 
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of the liquid sampling valve, it was established from the 
chromatograms that 2 - 2.5 microlitres of liquid sample was 
optimum for a 6 foot POROPAK Q column for the methanol-water 


system. 


However, acceptable measurement reproducibility was not 
obtained despite much effort directed to defining a G.C. com- 
puter job to analyse the peaks for the particular column 
Operating condition. The gas chromatogram was analysed using 
an IBM peak analysis package, which defined peak area calcula- 
tions based on first and second derivatives of the analog input 
Signal. A computer program was written to enable the user to 
obtain a graphical presentation of the raw input data and the 
first and second order derivatives. These results for a low 
concentration mixture are shown in Figure 5.5 and revealed 
the need for an amplified signal. The maximum output signal 
from the detector filaments was approximately 50 millivolts, 
and without amplification, excessive noise was introduced as 
the signal was transmitted from the G.C. to the computer, and 
in digitizing the small analog signals. The difficulty en- 
countered in defining exactly the beginning and end of the 
methanol peak, from the first derivative values shown in 
Figure 5.5, is obvious. This uncertainty, causing a varying 
methanol peak elution time was taken to be the source of the 


poor measurement reproducibility. 


To overcome the low signal readings, an amplifier with 
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en amplification of 30 was installed immediately after the 
detector filaments. The effects can be seen in Pigure 5.6. 
For a sample similar to the one producing the results in 
Figure 5.5, the derivative of the amplified signal has much 
less noise, decreased by a factor of 30 or more, and clearly 
defining the beginning and end of the methanol peak. With 

the amplifier installed, the parameters for the G.C. job de- 
finition for the computer analysis program became less cri- 
tical, and reproducibility was increased to a very satisfactory 
level. In essence, the amplifier reduced the effect of noise 
in the input signal, electrical noise due to interference bet- 
ween the computer and the gas chromatograph and noise due to 
poor digitizing at low signal levels by the analog to digital 


converter. 


As mentioned previously, three sample valves of diffe- 
rent design were used for taking a liquid sample and intro- 
ducing it in vapor form to the separation column. These were 
the Beckman slider valve, a Valco Hewlett Packard Option 513 
automatic sampling valve, and an Applied Automation Model Ix 


valve. 


The two main requirements of the valve were adequate 
flow of sample liquid through the valve and instantaneous 
flashing of the sample. Of course reliable, steady operation 
at high temperature over a period of time was desirable as 


well. The Beckman valve, with sample volume and a larger 
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bypass drilled into a sliding teflon block, was originally 
installed in a Beckman Series 4 Industrial Chromatograph, in 
the same oven and thus at the same temperature as the separa- 
tion column. This configuration had been used for distilla- 
tion column studies at the University of Alberta until 1975. 
Flow and operating reliability were adequate, but flashing 

of the liquid sample was slow, resulting in less than optimum 
separation. Adjustment of column conditions was difficult, 
as a high flash temperature but low column temperature are 
needed for good analyses. Also, temperature control was a 
trial and error procedure, with only a crude heater to adjust. 
At the commencement of this work, a new chromatograph with a 
separate valve oven and new sample valve were installed. The 
sample valve, manufactured by Valco, consisted of a conical 
teflon block with sample volumes cut into the sides, rotating 
in a metal sleeve with ports. The heat transfer and flashing 
properties of this very compact valve were much improved, as 
one face of the sample volume was metal. Flow through the 
valve, however, was not adequate, due to the complexity of 
the flow circuit. Breakin procedures were difficult, and the 
teflon block would scorch onto the metal sleeve with con- 
tinued use. Similar problems with a replacement valve of the 
same basic design led to the decision to search for another 


type of valve. 


The valve presently installed is an Applied Automation 
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Model IX. It is a unique design, in which pneumatically 
actuated small plungers move small fractions of an inch to 

block or open flow passages of carrier gas or liquid sample 
through the sample valve or bypass. The sample volume is 
encased in metal on 5 sides. Flow rates through the valve 

are higher than both -other valves and it has been used 
successfully during most of this work. A small teflon diaphragm 
must be replaced periodically due to deterioration with tem- 


perature. 


The sample valve was actuated under control of the IBM 
gas chromatograph package, which also performed the peak 
identification and area calculations. The composition in 
weight percent methanol, calculated from the areas determined 
by the IBM package, was then written into a DDC measurement 
Ioop by an online user written program. This program, desig- 
nated as DWL44 (see Appendix B), was queued by the IBM package 
at the end of the job cycle. Thus, the composition was 


available to the bottoms composition control loop. 


See. 3, Holdup 


An attempt was made to determine the mass holdup of 
liquid on each tray of the distillation column and in the 
reboiler and condenser for comparisons with values predicted by 
Common correlations. Furthermore, an attempt was also made to fin 


unique tray holdup correlation for this distillation column, 
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based on liquid and vapor flow rates predicted by the model. 


Holdup could not be measured but was estimated as follows. 
With one tray as a standard (feed tray), estimates of the 
relative quantities of liquid on the trays were made by view- 
ing the froth nature’ and froth heights through the clear glass 
Walls Of the column. The total weight of solution in the 
column was determined by operating the column at steady state 
and abruptly then stopping all flows, and draining the liquid 
from the column for weighing. The condenser and reboiler 
holdups were determined with the column non-operational. 
Since the levels of the condenser and reboiler are controlled, 
it was assumed that their holdups did not change so by sub- 
tracting the condenser and reboiler holdup from the total 
column holdup, the cumulative tray holdup could be calculated. 
Individual tray holdup was calculated by taking the relative 
volume estimates, applying a density correction, and multi- 
plying this|tray ratio by the cumulative tray holdup. This 
procedure was followed for the standard operating conditions, 
and for feed flows of ten and twenty percent above and below 


the normal flow rate. 


The results of the estimated holdups are presented in 
tabular form in Appendix F (Table F.1) and in graphical form 
in Figure 5.7, for all operating conditions tested. Also 


plotted on Figure 5.7 are holdups based on published holdup 
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correlations, using liquid and vapor flow rates from indivi- 
dual trays, predicted by the model described in Chapter 3. 
Both experimental holdups (for a variety of operating condi- 
tions) and holdups derived from correlations predict a similar 
trend of holdup decrease with an increase in tray number. 
However, holdup correlations predict values that are larger 
than measured ones, and in the case of the Francis Wier pre- 
diction, larger by up to a factor of two. A complete summary 
of estimated and measured tray holdups is presented in Appen- 


ex: F 


Attempts were made to correlate estimated holdup values 
with both tray liquid flow rate (predicted by the model des- 
cribed in Chapter 3) and feed flow rate, using the plots shown 
in Appendix F. No definite relationships were observed in 
either case, although generally, an increase in feed rate 
appears to result in an increase in holdup on all trays. The 
excessive scatter of data is attributed to the estimation 


technique, and may be observed in Figure 5.8. 
The published holdup correlations used were as follows: 


1) AICHE’? method of efficiency prediction for bubble cap 


Erays 


tam 1.65 + 0.19w - 0.65F + 0.020L ares 
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where ie PnChesrOfwliquid onthe tray 


w = weir height (inches) 
O35 
F = V Pa 
Me = Vapor tlow fLrom tray (ft/sec) 
®g = Vapor density (lbm/£t?) 
Pee = togu.d tlow (gpm/it) of flow width) 


Volume holdup was converted to a mass by applying a 


density. 


- pate 
1i) The Francis Wier Sore for flow over a rectangular 


wier: 
spies: 
C= .415:(0 = 0.2h AS v2g GS 2) 
: | 
where q = volumetric flow over the weir (ft /sec) 
L = length of the wier crest (ft) 
AS = height of liquid above wier (ft) 


; 2 
fecal vacceleration due to. gravity (ft/sec) 


This equation was used to calculate a weir head for a 
given flow through the weir. The assumption that the term 
ec. 2h © was negligible next to the wier crest length L simpli- 
fied the equation. The entire liquid volume on the tray was 
assumed to be clear (negligible froth volume), and volume and 


weight holdup was calculated on that basis. 


iii) Foss and Gerste?-holdup equation 


oN: 


hy = 0.24 + 0.725h, - 0.29h, v_ ae Pee Agha) ziteN5e3) 
where ne = outlet wier height (in) 
he = liquid head on tray (in) 
em = gas velocity (ft/sec) 
Pg = gas density (1bm/£t?) 
q = liquid flow rate (ee Jae) 
Z = average liquid flow width (ft) 


Sample calculations for all cases were made in, or the re- 
Sults converted ito,SI units, and are also presented in 


Appendix C. 


Holdups predicted by the published correlations are all 
higher in magnitude than those measured. This may be explained 
by the fact that bubble cap trays, designed and constructed 
by Svrcek 7° fora the, original distillation column control 
study, may not be of a standard design. Also, the column is 
a pilot plant column, smaller in scale than industrial columns, 
for which the published correlations were derived. For this 
reason, these published holdup correlations were not used 
for predicting holdups during the computer simulations of the 


distillation column, as presented in Chapter 4. 


2.2 RESULTS 


Considerable time was spent in determining an appropriate 


set of operating conditions. It was found that a tradeoff 
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between sensitivity to disturbances and bottom composition 
measurement accuracy was to occur. It was decided that pro- 
duct compositions should be between 4 and 5 percent in the 
bottoms and between 96 and 97 percent methanol at the ejey 
The composition range was similar to that used by other 
workers, but the column had been operated at lower bottoms 
composition to achieve good distillate sensitivity to upsets. 
For this reason, certain disturbances will exhibit lower 
gains than open loop responses from previous tests by other 
workers. A list of sample initial operating conditions can 


ee cound in Figure 5.9. 


As previously mentioned, a gas chromatograph was used to 
measure the bottom product composition. This is a discrete 
measurement device, and inherently has a dead time of 4 
minutes (in the analysis program used in this study). As 
will be noted in the plots of bottom composition responses, 
there appears to be a time difference of from 8 to 12 minutes 
between the simulated and experimental trajectories. This 
may be explained with the aid of Figure 5.10. The G.C. pro- 
gram of the IBM 1800 computer was used to collect gas chroma- 
tograph analyses. Upon completing an analysis, a G.C. queued 
service program (DWL44, Appendix B) would write the analysis 
result into a measurement word =f DDC copy recora.” ~ve= 
pending upon when the nese aenkl was taken with respect to 


the start of the disturbance, the delay in getting a measured 
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FIGURE 5.10 Bottoms Composition Analysis Collection 
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composition into the DDC loop record could be from 4 to 8 
minutes. Another data accumulation loop with a pole time of 
256 seconds would take the measurement from the measurement 
loop, and write it to a measurement history file. If this 
loop were not in phase with the gas chromatograph analysis, 
an additional delay of 256 seconds could result. If the 
measurement were filtered, another delay could be 

Bncurred. Together, this could constitute a total of 16 
minutes delay. However, all response lags fall within the 


4-12 minute delay, as filtering was not used. 


There also exists a slight process delay due to the G.C. 
Gycle loop (section 5.1.3, Appendix E), but this delay was 


considered to be small compared to other G.C. delays. 


Results are presented for positive and negative 20 per- 
cent step disturbances in reflux and feed rates, and 10 per- 
Bene Steps in steam rate to the reboiler. In Figures 5.11 
to 5.16 the experimental response in each case is accompanied 
by a simulated response. The simulation results presented here 
are for the 04 series (Table 4.1), with the assumptions of 
constant mass holdup, variable heat loss and variable effi- 
ciency, with the variation according to changes in top compo- 
sition. This simulation series was chosen because it was 
found to adequately represent the experimental response. To 
ascertain the agreement that would exist between predicted 


responses based on different assumptions and the experimental 
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FIGURE 5.11 +20% Feed flow; Simulation vs Experimental 
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FIGURE 5.12 -20% Feed flow; Simulation vs Experimental 
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FIGURE 5.13 +10% Steam flow; Simulation vs Experimental 
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FIGURE 5.14 -10% Steam flow; Simulation vs Experimental 
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5.15 +20% Reflux flow; Simulation vs Experimental 
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FIGURE 5.16 -20% Reflux flow; Simulation vs Experimental 
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results, the reader should refer to Chapter 4. Some of the 
initial experimental steady states vary Siigiecly Crom run to 
run, varying by +.05 percent methanol, which is close to 
instrument resolution. These runs were taken over a period 

of weeks and although every attempt was made to keep the column 
Operation constant, such slight differences appeared and were 
considered normal. In simulations, no adjustments were made 


for this minor difference in composition. 


ao DISCUSSION 


Comparisons between experimental and simulated responses 
to disturbances may be made in two distinct areas to assess 
the validity of the model used for the predictions. The 
Pirst is prediction of terminal compositions at steady state 
of various operating conditions. The second area for compari- 
son is the column dynamics or rate of change of column com- 


positions due to changes in operating state. 


The experimental responses themselves are representa- 
tives of duplicate runs. Taken over a period of approximately 
two weeks, they are deviations from a set of reference steady 
state conditions decided upon with criteria of column sensi- 
tivity and composition measurement accuracy in mind. It may 
be noted that top composition response to a step increase in 
feed flow rate is almost undetectable. This was the result 


; ; ; ++ 
of a compromise between column sensitivity at different 
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operating conditions, and instrument Sens wtivacy. 


The steady state conditions resulting from the moderate 
disturbances in steam, feed and reflux rates from the reference 
steady state values agree well with the experimental behavior 
in all cases. The maximum discrepancies may be found in the 
feed disturbance runs, being of approximately 0.2 percent 
methanol (or 5% of range) in the bottom composition for both an 
increase and decrease in flow. The steady states here have 
been matched more closely by taking greater care in choosing 
heat loss and efficiency fitting parameters. In these cases, 
as in all, heat loss was taken as uniform for all trays, 
whereas tray efficiency was considered uniform throughout the 
Stripping section and again throughout the rectifying sec- 
tion, but both heat loss and efficiency were allowed to vary 
with terminal composition to match the experimental steady 
seater.) Although, at was not possible to.determine, if. the 
efficiency was representative of the experimental value, the 
heat loss chosen to match the compositions at steady state 
compared well to the heat loss from the physical column as 
calculated by the energy balance program BALANC (see Appendix 
B). The tray efficiency, also used as a fitting parameter, 
was found to be less effective in this regard than heat loss 
and could not be determined by measurement but only calculated 


from column operating data. However, the’ fitted efficiency value 


used in the simulations ranged between reasonable limits as 
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can be seen from Table 4.3. 


The transient response dynamics agree very well in all 
cases for the bottom product and top product composition 
transients involving reflux rate increases and steam rate 
decreases. However, the feed disturbances, and other distur- 
bances causing a decrease in methanol concentrations in the 


distillate require further analysis. 


Unique to the feed disturbances is an initial resvonse 
in a direction opposite to that expected and to the eventual 
trend. This is especially prominent in the negative feed 
disturbances, but always only in the distillate response. 
Although no simulations were run testing the effect of varying 
efficiency on the stages to vary intermediate compositions, 


it is expected that the explanation may be found there. 


For simplicity, efficiencies were taken as uniform from 
ima tOutray at some value in the stripping section), and at 
another throughout the rectifying section. This simplifica- 
tion was used by Stainthorp and Searson* in a transrer siine. 
tion model study. In this work, a trial and error procedure 
was used to match the initial steady state. Once this steady 
state had been matched, these initial efficiencies were used 
in all simulations. As can be seen from a sample steady 
state presented in Figure 5.17 the resulting composition pro- 


file had a feed tray composition of 51.0 percent methanol. 
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Thus, a methanol balance about the feed tray shows that the 
input flow with a composition of only 50.0 percent methanol 
tends to lower the composition on the feed tray. When this 
input term is reduced in magnitude, the effect would be to 
let the feed tray composition rise momentarily, sending a 
Besse UD.and down the column. It is postulated ‘that-this is 
the effect that manifests itself in the response shown by the 
distillate, the same effect being damped out by the large 
reservoir in the reboiler. This inverse response phenomenon 
was also observed by Luyben and Vinante*> on avpiloteplant 
column, but was unexplained, as the study dealt with trans- 
fer function models for terminal decoupling. By the same 
token, a negative response could be predicted for a step in- 
crease in feed, under the same initial conditions, as is also 
observed. These inverse deviations are very slight, only 

0.1 percent methanol at the feed tray, but it can be seen 
from the simulation results that this pulse magnitude is 


maintained through the rectifying section. 


To check this explanation further, it would be necessary 
to use different values of efficiencies that would predict 
identical terminal compositions, but different intermediate 
compositions at steady state, more specifically, a lower com- 


position on the feed tray. 


Also showing moderate deviations in dynamics from the 


experimental responses, and somewhat peculiar curves in the 
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response, are disturbances resulting in a decrease in dis- 
tillate composition. Duplication of the response with an 
alternate integration method and integration time step eli- 
minated the stability of the integration as the fault... A 
behaviour such as exhibited by the distillate composition 
under a negative feed flow disturbance would indicate a series 
of first order systems of varying time constants between the 
geea trayoand the distillate tray... In a matrix system, it 
would indicate varying eigenvalues in the "rectifying sec- 
tion" of the matrix. These eigenvalues are completely a 
function of liquid and vapor flow rates and compositions, and 


would change as compositions change in a nonlinear model. 


The peculiarity of some of the responses in question may 
also be attributed to the choice of efficiency. As previously 
mentioned, efficiencies were arbitrarily taken as uniform 
throughout the stripping and rectifying sections. Certainly 
this forces the simulated tray compositions and flows to 
assume an arbitrary profile; as vapor and liquid flows, liquid 
heads and compositions vary from tray to tray, so will actual 
efficiencies. Each tray may be regarded as a buffer zone at 
steady state composition, with circulation between buffer 
zones on either side. There are composition gradients from 
the reboiler to the condenser and there are circulation rate 
gradients, with a discontinuity at the feed plate, as the 


feed is injected there. A buffer zone with a composition of 
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concentration gradient out of line from the natural gradient 
would retard or propogate a disturbance in an unnatural way, 
as its eigenvalue changed in an unnatural way. A changing 
efficiency or heat loss causing a variation in an eigenvalue, 
would cause a change in the rate of change in the column. 

For example, the heat loss of tray 6 physically changes ac- 
cording to tray 6 temperature, but in the simulation, heat 
toss is varied with top composition. This may cause the tray 
6 liquid flow to react more slowly than in the physical 
column, and cause the simulated composition to also react 


more slowly. 


All simulations presented with experimental data here 
are for the case of constant mass holdup. Also, efficiencies 
and heat loss were taken as functions of distillate composi- 
tion. For disturbances resulting in a decrease in distillate 
compositions, the dynamics appear to be faster than shown by 
the actual column. There may be several reasons for this. 

In the actual column, holdups vary with changing conditions, 
mSeholdups are a determining factor in tray dynamics, changing 
Conditions affect column dynamics; the effects of varying 
holdups are not presented here, but in Chapter 4 of this 

work. Also, tray heat losses were varied as functions of 
distillate composition, but in actuality, are functions of 


tray temperature. In every type of disturbance, the heat 


loss on a tray would naturally rise as tray temperature rose, 
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reducing the magnitude of the disturbance to the next tray, 
ma slOWing it downs Again, it was only for simplification 
purposes that this assumption was made. The model solution 
scheme would readily accept tray heat loss as a function of 
tray temperature, and only through further testing can the 
effects of that simplification be established. The same holds 
eeue for tray efficiency, as efficiency was taken as a func- 


tion of terminal composition also. 


It is interesting to note that only in responses in- 
volving a decrease in distillate composition, are these irre- 
gularities exhibited. The large holdup of the reboiler, five 
times that of an average tray, effectively damps out the 
Beignt deviation for all types of disturbances in the bottom 
meomoositton “transient. In the distillate, disturbances causing 
a decrease in distillation composition, because of basic non- 
Minearities of the distillation column are typically more 
severe than other kinds. For increases in distillate compo- 
Sition then, stronger resistance to changes, or stronger 


damping, eliminates the effects of peculiar dynamics. 


To test these explanations further, testing with alter- 
nate efficiency profiles and heat loss correlations must be 
done. The assumptions of uniform efficiency and heat loss 
may then be shown to be too restrictive. However, this 
testing without physical measurements to reference is a trial 


and error procedure and thus time consuming and expensive. 
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The simplifications chosen, however, do not limit the gene- 
rality of the solution scheme of the model, as any efficiency 


or heat loss profile could have arbitrarily been chosen. 


Despite these simplifications, the model predictions com- 
pare very well with the experimental data at every operating 
steady state. In all cases, the bottoms compositions tra- 
jectories simulated are very close to those obtained experi- 
mentally. Top composition responses for steam decrease and 
reflux increase disturbances predicted by the model match 
very well with the experimental responses, and for other dis- 
turbances within a reasonable error limit. The solution 
scheme of this model was thus considered valid, and was then 
used in comparison studies of models of various degrees of 


Simplification. 
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CHAPTER 6 


CONCLUSIONS AND RECOMMENDATIONS 


6.1 CONCLUSIONS 


A nonlinear solution scheme to model a binary distil- 
lation column and to include the effects of variable mass 
holdup in the model has been developed and tested for distur- 
bances in feed, reflux and steam flow rate. Comparisons to 
actual experimental responses of a pilot plant column have 


shown favorable agreement. 


Simulations were performed with the rigorous model altered 
using several different common simplifying assumptions. The 
predicted behavior was compared to the prediction of the 
rigorous model and the effects of the various simplifications 


were made evident. 


To obtain good agreement with experimental responses, 
hE was found that fitting of simulated responses was an inte— 
Gral part of the modelling effort. Tray heat loss and tray 
efficiency were used effectively as fitting parameters, the 
column being more sensitive to the heat loss parameter. It 
was also found that overall column dynamics were strongly 
influenced by the dynamics of the fitting parameters. It 
appeared that the dynamics of the reboiler composition were 


largely determined by the reboiler time constant, other time 


constants being relatively less significant. 
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The model was found to be stable in the solution of a 
nonlinear equation set including the holdup correlation 
equations. The hydraulic dynamics were found to be the sta- 
bility limiting element, composition dynamics being much 
Slower. Three types of holdup correlation equations were 
used in variable mass holdup modelling, and it was found that 
overall column dynamics could be altered in both directions, 


depending upon the holdup correlation parameters chosen. 


Finally, the solution scheme developed in Chapter 3 has 
been shown to be creditable, its solution satisfactorily represen- 
ting actual physical data for disturbances studied, and is 


easily adaptable to simplification if desired. 


6.2 RECOMMENDATIONS 


It is recommended that effort towards further distilla- 


tion column modelling be directed toward the following: 


1. Elimination of the fitting parameters. It would be 
useful to have heat loss and efficiency data in 
functional form, relating these parameters to ope- 
rating variables. 

2. To further experiment with different holdup corre- 
lation parameters, and through simulations attempt 


to further define their effects upon simulated column 


dynamics. 


2k 


In addition, a problem exists in the data gathering of the 
bottoms product analyses, which should be remedied for further 
distillation modelling studies. Gas chromatograph analyses 
must be properly phased with disturbance times and DDC accumu- 
lation loop pole times, to reduce the time delays observed in 


eaas work. 
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NOMENCLATURE 

mo, c, a - parameters, unique to particular equation 

A - system matrix, defined by composition balance 
equations 

A - intercept of enthalpy - steam temperature 
relationship 

B - slope of enthalpy - steam temperature 
relationship 

D - distillate flow (g/sec) 

DT - integration time interval (sec) 

E - Murphree Vapor Efficiency 

ply = vec cOhvused in solution of vapor and) liquid 
flows 

A = vapor enthalpy of vapor leaving stage n (J /gm) 

Do sv liquid enthalpy of liquid leaving stagein 
(kJ/gm) 

HSI - reboiler inlet steam enthalpy (kg/gm) 

HSO - reboiler outlet steam enthalpy (kg/gm) 

i - identity matrix 

K =~ ratio iofi vapor to liquid composition on each 

: Gray 

Ln - liguid flow rate leaving stage n (gm/sec) 

5 - parameter (heat loss or efficiency) 

Q - heat input, heat loss (kg/sec) 

R - reflux flow rate (gm/sec) 

S - Laplace domain parameter 

S ~ sidestream flow rate (g/sec) 

S -~ steam flow rate to the reboiler (g/sec) 


4 a 


‘id 


Nn 
O 


Subscripts 
o 


i 


Superscripts 
K 


* 


Nomenclature 


AHLO 


aA) 


condensate flow rate from the reboiler (g/sec) 
temperature (°C) 
time (sec) 


matrix used in solution of liquid and vapor 
flow rates; has enthalpy and holdup terms 


heat transfer coefficient in reboiler 
(kJ /sgec°c) 


vapor flow rate leaving stage n (gm/sec) 
liquid holdup on stage n (gms) 
liquid composition (% methanol) 


vapor composition (% methanol) 


condenser 

initial 

final 

reference condition 
stage number 


reboiler 


denotes variable value at the Kth time instant 


equilibrium value 


intercept for steam enthalpy correlation 
(kJ/gm) HS = (BHLO * T + AHLO) 
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BHLO -~- Slope for steam enthalpy correlation (kJ/gm°C) 

eal - liquid coefficient in variable holdup corre- 
lation 

(ey - vapor coefficient in variable holdup correlation 

C3 - intercept in variable holdup correlation 

D - (in ENTH) derivative of liquid enthalpy with 


Becnmect to. Liquid composition 
- (in MAIN) disturbance vector values (gm/sec) 


DHL = derivative of liquid enthalpy with respect to 
time 

DT - integration time interval (sec) 

DWT - derivative of holdup with respect to time 

DXT - derivative of liquid composition with respect 
to time 

DTL - derivative of liquid temperature with respect 
to time 

D2 - derivative of tray temperature with respect to 


licusds COmpOsiELOn 


EE - efficiency at the initial steady state 

E - efficiency at final steady state 

EF “Clr reLlengyvsat), current. Cime 

F - feed flow rate vector (gm) 

FL - initial feed flow rate (gm) 

BT - (in INTGR) work vector - ratio of vapor to 
liquid composition, used to obtain tri diagonal 
matrix 


- (in VFLOW) work vector used in the solution 
of flow rates 


HF - feed enthalpy (kJ/kg) 
HLT - liquid enthalpy (kJ/kg) 


HVT - vapor enthalpy (kJ/kg) 


Ye 


fade ‘aL OE 2 tinh oP kB: am " 


: | 7 r 
aay ac kaw +): tequ 


tae. Saiw QC TeORes nian 


cine) aeneinh whe. leicht uted 
vi rabye es teat ian Bie wie 74Ne eRe < “Oi 
oxwajet's “tieiN fois secgian i sl Ms a wd OF ; 


rf 7 
‘ ~ 
As woes sr aiiad te Sy. sev ca 


by 
63 spegney see - aie paseesaa ie 


i A baw base * jana nt 


+ 


ISTRN 
IES T 
JM,JN 
KK 


KT 


BT 
LS 
MTT 


NTT 


QLP 


ROG 


liquid enthalpy data vector (kJ/kg) 
vapor enthalpy data vector (kJ/kg) 
steam enthalpy (kJ/kg) 

vector of disturbance times (min) 
disturbance indicator 

printing counter, PRINT index 
printing counter, PRINT2 index 


vector of possible disturbances 


ASCII characters indicating disturbance type 
variables specifying print times (sec) 


A eyolewindicator 


gL Ps0 


Imcica tox orcoperati ons joccurring (on tne earse 


program cycle (eg. read data) 
liquid flow rates (gm) 


initial feed flow (gm/sec) 


number of trays plus reboiler plus condenser 


numberof trays plus reboiler 
number of trays 

simulated run time (min) 

heat lossrat vinitralesteady tstate 


reboiler heatwinput (Cj7sec) 


heat loss at final steady state (j/sec) 


heat loss at current state (j/sec) 


initial reflux flow (g/sec) 


ratio of feed to feed tray composition (to 
obtain a tri diagonal matrix representation) 


vapor density (g/l) 


ROL 


TA, TSU 
TLS 
TR 
UA 
VT 
WT 
XF 
ALT 


XIF 


XIN 
AT 

XST 
XXT 
XBT 
aay 

YST 


YS 


liquid density (g/l) 

Steam flow rate (g/sec) 

Steam flow rate at initial steady state 
steam temperature (°C) 

time (min) 

tray™temperature vector (°C) 

base matrix in integration scheme 
Matrix in integration scheme 
temperature-composition data vector (°C) 
reboiler liquid temperature (°C) 

heat transfer coefficient (j/sec°F) 
vapor flow vector (g/sec) 

holdup vector (gm) 

feed composition vector 


initial feed composition 


final composition interpolation value (for 


varying heat loss, efficiency) 

initial composition interpolation value 
liguid composition vector 

liquid composition data vector 

prenting butter for top composition 
printing buffer for bottom composition 
vapor composition vector 

vapor composition - enthalpy data vector 


vapor composition data vector 
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APPENDIX A 


THE METHANOL WATER SYSTEM 


TABLE A.1l EQUILIBRIUM, TEMPERATURE DATA 


X Ng T, G) T. C) 
0.0000 0.0000 EOO020 POOR” 
Die 3 50 Os2153 9373 O97.6 
0.0689 0.3469 9520 O97 re 
0.2020 0.4372 91.4 Shee 
pet 339 Posy claps) Soa Ise 
0.650 O56 09 Salta) SMe: 
02026 0.6143 = - 
052389 Dede ayers 84.3 96.4 
One739 0.6870 = = 
U0 078 O 7100 Sse 94.9 
0.3405 Oia Sa vs = 
O10 2 0. 7500 1956 93.2 
0.4028 Orahio 32 = = 
0.4325 8 er er AS! TIS 956 
0.4612 On7 927 = = 
0.4892 0.8043 = 5 
Oo L62 OAs wlesys] = a 
0.5246 Oris 271 ieee) 88.8 
OW5926 0.8449 = a 
0.6400 0.8624 one S47 
D293 0.8934 fi.3 80.8 
0.8508 0.9224 69.4 76.6 
Oro 167 0.9504 ei re, Zee 
0.9412 029759 eieye Js 685.3 
Ces 2 ae ate Oio.6 oa. 65.4 66.4 
1.0000 1.0000 64.7 64.7 
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TABLE A.2 


oaoKt Kx 


X (wt) 


0.00 
0.05 
U0 
OFet> 
O20 
Oix25 
0230 
On0.0 
0.40 
0.45 
0.90 
Oya5)5 
0.60 
0.65 
Le, 
Oar fs, 
0260 
0.20 
00 


liquid phase 
vapour phase 
liquid phase 
vapour phase 
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we sS 
-4370 
yI342 
-6140 
~06056 
e710 
-(3G2 
HOD 2 
7874 
. 8087 
wooo) 
-8476 
moO Og 
Seo oy 
2708 
pees 
O26 
0000 


oOo Oro Oro (Or OO OS aa Os Ore 


methanol composition 
methanol composition 


enthalpy 
enthalpy 


(kj/kg) 
(kj/kg) 


EQUILIBRIUM, ENTHALPY DATA 
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PHYSICAL PROPERTY CORRELATIONS 


The following correlations were converted for use with 
variables in the SI unit system, from correlations for use 


with the engineering system of units. 


a) density of liquid water 


Pomme) =) 005s =22414 x 107¢¢ ae ee Loar (Agia) 
a0 eee Te <3 120°C 
b) density of liquid methanol 
B(gm/ml)-='0.808 <071740x 1077 - 1.98 x 10° °n* (A.2) 
On Ceca) = 1205.0 
c) density of methanol-water mixture 
s a =3 - 
Be Cay SO 20 Sell 0 1arsdeSdk10 °T"-1.52x101 x-e5) doe ee 
a0 7CheeT < 120% (A.3) 
QO wt% MeOh < X < 100 wt% MeOH 
d) density of saturated methanol-water solutions 
= -6.,2 
ecm mi 0" 960) = 11 36. LO khan hed Gux 1 Oma am eee 


O wt3 MeOH < aS 100 wt% MeOH 


e) density of steam 


ee 


>m = 4,299x10° 'T* + 


2 ag 10 LO 


9(gm/om?) = -1.018x107 
Sy Sep gh ao (A.5) 
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f) heat of vaporization of water 
H(kd/kg) = 2487.4 - 1.826 x 107-7 - 4.24 x 107¢r? (a6) 


6620) <\T¥s, 150°C 


g) heat of vaporization of methanol 


H(kJ/kg) = 651,.- 1.80T - 1.156 x 10 “fT (AT 


h) heat capacity of steam 


Breen Tome LO ioe ee ee 


a acre 7428) se Von ete oe om ee 


= eal/gm°c 


NARA Le mak @ 
5 /kg 


i) heat capacity of liquid methanol 


¢, cal/gm°C = .557 + 1.87 x 100 1 9,01 eoloee ns 
c, kI/kg °C = 2.33 + 7,82 x LOG2T ATT Cie 1 ee 
j) heat capacity of vapour methanol 
is ene 
c, cal/gm°C = .3206 + 6.19 x 10 ones) 45 ex OPE 


We oO eae em 


ewes Ou On st 8G oe OS 


oP LOT Sa PS ae amd @ 
> /kg 


2.5 x07) 9 n? (A. 10) 


k) saturated liquid temperatures of methanol-water mixtures 


-2.2 
mee Oy =,99.6 — .861X + 1.02 x 10 x 


oh ree saa oe (xan 


O wt% MeOH < X < 100 wt% MeOH 
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1) saturated vapor temperature of methanol-water mixtures 


Tec = 99 op e-h p25 xo ng tx 2059 yR eto x2 


Br OSARe LOX? (A.12) 


O wt% MeOH < X < 100 wt% MeOH 


m) enthalpy of saturated methanol-water vapor 
Pees KG) = 2655.3 —- 12.29xX (A.13) 


0 wt% MeOH SUG EES 100 wt% MeOH 


n) enthalpy of saturated methanol-water liquid 
-2.2 


cd ZG) = 41357 = 5.75X + 3.15 x 10 “x (A.14) 
O wt < X < 100 wt 
©) vapor pressure of liquid methanol 
oo Hoymae 07 30° —TLe sa t7 (T{ °C) 250) (A.15) 
Cee Leo ZU € 
Pp) vapor pressure of liquid water 
aectig po Hg) = 7.9668 - 1668.2/(T+220) (A.16) 
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APPENDIX B 


ON-LINE AND SIMULATION SOFTWARE 


Computer programs used in this work were either used in 
gathering transient data for open loop runs, steady state 
information gathering and reporting at the beginning and end 
of the experimental runs, or were used in the distillation 


samulations. 


Data accumulation and control of all pertinent variables 
on the distillation column was handled by an on-line IBM 1800 
computer with a DDC package. All control or accumulation 
loops were accessible through a teletype near the pilot plant 
column. Other loops provided the capability of storing his- 
feorcal information. A program RETRV could access the his= 
Sorical data, type it on a terminal or punch it onto cards 


for later plotting. 


The programs DASS, DATAC and BALNC were used to measure 
and report column conditions at steady state. These programs 
were run at the beginning and end of experimental runs to 
check calibrations of flow and composition measurements by 
mass, composition and energy balances. These programs were 
modified from the originals used by other workers at the 


University of Alberta in that a conversion was made from 


Engineering units to the new SI system efsunits. in ‘alle caica- 


lations and report formats. 
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The gas chromatograph system involved the use of several 
eine “computer programs; G.C. monitoring and peak analyses 
were implemented by the IBM G.c. program package available 
on the IBM 1800 computer used to control the pilot plant 
column. The package enabled automation of the bottoms pro- 
duct sampling system, and in conjunction with a program DWL44, 
used peak areas to calculate methanol composition, and pro- 
vided output on a teletype with each analysis. A special 
DDC data accumulation loop was updated with each analysis. 
The GCRAW group of programs, specially supplied by the DACS 
Semele start aided invG.C. job definition; the programs were 
Used to scan G.C. output of a typical methanol-water chroma- 
togram and plot the raw data along with calculated first and 
second derivatives. The plots were then used to choose 
realistic G.C. job parameters. Documentation of the GCRAW 


programs is available from DACS centre staff. 


All distillation column simulations were run on an IBM 
360 and later an AMDAHL 470/V6 operated by the Computing 
Services department at the University of Alberta. Simulation 


programs were kept on disk files and were run from a remote 


terminal. 


The source module for the main modelling program was 
kept under US010.S, and a listing is provided in this Appendix. 


Modification to implement simplifying assumptions were made 


to this program, but copies of modified programs were not kept. 


PY Ie 


: a H0R 
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Data required by the simulation were kept in three 
files, one for each type of disturbance examined. The names 
memo t cd, US020.d, USO30.d for reflux flow, steam flow 


and feed flow disturbances respectively. 


Calculated output of the simulation program was kept in 
several different scratch files with varying formats. When 


desired, the files were output to cards, for later plotting. 


The program EPLOT was used to plot transient column res- 
ponses, both experimental and simulated. Several variables 
could be plotted on one set of coordinates, for comparison 


purposes. 
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APPENDIX C 


DERIVATIONS 


The solution schemes described in Chapter 3 had the ob- 
jective of solving the system of 40 equations, or 4 systems 


of 10 equations, shown below, where "n" denotes the stage. 


a) composition balances 


dx 1 
dt a wr 'Tn+1*n41 ry fs Ea ee i Vee 
qwrt. 
a eee erial 
era ] ( ) 


b) total balances 


LOOT Ge 


Aree So KC Jz) 
ar n+l a “A Ua 


c) enthalpy balances 


dh 1 
ae = 2 Pan 
eke WT (he ns En james 
dwt. es 
a hy ar aat QLP] . 


@) holdup correlations 


= + V GP 
WT a_L b n 


These equations are not rigorous, but specific to this 
work; the limiting assumptions used in this model have already 


been applied. These equations apply to the general equili- 


brium stage. Terms must be added or deleted for special 


usZ 


Stages such as feed tray, reboiler and condenser. The colum 
meramnecers x,y,L,V and WT are all functions of time t. The 
following solution schemes attempt to solve this system of 


40 equations, although by two different methods. 


eee SUBSTITUTION METHOD 


In both the substitution method and matrix inversion 
methods, the coefficients of the composition balances were 
determined from the remaining three systems of equations. It 
was the composition balances that were directly integrated 
to obtain the trajectories of terminal and intermediate column 


compositions. 


The energy and total balances were solved together to 
obtain preliminary liquid and vapor flow rates. The co- 
efficients of these two sets of equations are enthalpies 
(directly obtainable from current compositions), holdups and 
holdup derivatives obtained from the previous time step. 
Holdups and holdup derivatives were recalculated in a correc=- 
tion step, where the preliminary flow rates were substituted 
into the holdup correlation equations (C.A4)om ethe new yhold= 
and flow rates were then used as coeffi- 


ups, accumulations 


cients in the next integration step. 


in solving for the 20 liquid and vapor flow rates, 20 


equations were used. Specified variables were reflux flow, 
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feed flow, reflux and feed enthalpies and heat input to the 
reboiler. The method for accomplishing this is outlined be- 
low, and mainly involves consecutive substitutions of expres- 
sions of variables in terms of each other, until an expression 


for one variable (Vg) in terms of known quantities is obtained. 


The systems of equations used to solve for liquid and 
vapor flow rates are Equations (C.2) and (C.3). The deriva- 
tive terms are directly obtainable from information from the 


previous time step, where K denotes the time step. 


K 
yonn = (wrk - wrk); 
at Dr n n 
ank see 
ee KC on Th 
at reat 


The coefficient Ch is the slope of the enthalpy-composition 


relationship, and continually varies with composition. 


The liquid flows are eliminated from Equations (C.3) by 


substituting equations from (C.2). 


stage 9 


h., - Lghg - VoHy - QLPg = WTgDhy + hgDWT, 


Vorg 7 10 92g 9 


Soe 810 
substituting for Ly 


Z ) 
VgHe + Lighyg - Bg (Vg * Ho 9 9 gHg 
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or 


Vg(Hg - hg) + Ly lhy, 9 


= WT,Dhy 


Following the same procedure for stage 8 


V>(H, - hg) + Lo(hy - hg) - Vg(Hg - hg) - OLP, 
= WT,Dh, 
but 
Lg = Vg + Lig - Vg - DWT, 
substituting 
pe temem ot th, = ,) it Dy yg Ae cove: 
- QLP = WT,DH, + (hg - hg) DWT, 


One can continue this procedure for all stages, including 
terms for the feed flow rate and the reboiler heat duty in 
the equations for stages below the feed stage and.) in Gene se— 


boiler equation respectively. 


If all non-vapor flow terms are collected on the right 
hand side of the equations obtained, and all terms in vapor 


flow on the left hand side, a matrix representation can be 


formulated. Let V be a vector of vapor flows and let F be a 


vector of terms on the right hand sides of the resultant 
equations. 
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+ Bi Baa] - hh) + WT dhe 
ar P(A i A) 
and 
9 
tn = QLP. + (hoa - h) 2 DWT 
i=n 
+ Li 9 (Ai -h) + WT Dh, 


and where the matrix A has the form 


a <b 
Dy 5 Cc 
re me 
Db. 
A = 
“9 
biagce 
Dea, 


mhe elements in matrix A are: 
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All terms in the F vector defined above are known, and the F 


vector may be fully calculated. 


defined as all enthalpies are known. 


The kh Mateixeisealso fully 


All liquid flow terms 


(except reflux flow Li 9) have been eliminated from the system 


of equations represented by (C.5). 


With further substitution, 


an expression for the vapor flow from stage 9 may be found, 


in terms of known quantities. 


Peommequations (C.5))<andy (C.8) 


stage l 
fy = 
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stage 2 
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ite 
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V — 
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for stage 
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This expression for V, may be solved, as all terms on the 


9 
right hand side are in known enthalpies and derivatives. By 
back substitution, all other liquid and vapor flows may be 


calculated. 


The change of liquid and vapor flow rates is accompanied 
by a change in tray liquid holdups, and the holdup correla- 
tions (C.4) were used to calculate these new values. Rates 


of change of holdup were then calculated by Equation s(ey9) 
Sie x = (G59) 
DWT (wTs WT) /DT 


where the superscript * denotes the updated holdup and deriva- 
tives. Before completing the cycle by. takingsthetupdatedsco- 
efficients to the integration step, final corrected flow 

rates were calculated via the substitution method just des- 
Cribed. It was found that one such iteration was sufficient 
to accurately correct flows for changes in holdup, and vice 


versa. These updated values of holdup, rate of holdup change, 


and liquid and vapor flow were used in further calculations 


in the next time step. 
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With new holdups, flow rates and derivatives, the co- 
efficients in the composition balance equations had been 
determined. These equations were then integrated via the 


modified Euler or other integration scheme. 


In summary, with the substitution method, the energy 
balances were presolved together with the total balances for 
the general case, yielding expressions which were solved for 
liquid and vapor flows. Holdups were calculated from liquid 
and vapor flows, and the latter were again corrected for 
holdup changes. The new holdups, flow rates and derivatives 
were used in the integration step as coefficients, to complete 


the integration cycle. 


C.2 THE MATRIX INVERSION METHOD 


Unlike the substitution method, the systems of equations 
were formulated into a matrix format, and were solved by com- 
puter library routines, by matrix inversion and vector opera- 
tions. Other differences are evident in the choice of the 
systems of equations that are solved together. Ini ithe matrix 
inversion method, the energy balances (C.3) and holdup cor- 
relations (C.4) were solved together to calculate vapor and 
liquid flows, and the component balances (C.1) and total 


balances were solved by integration. The tray hoidup deriva- 


Pee tearm was eliminated from the energy balance by substli- 


tuting in the total balance (C.2). Equations (C.10) and 


Log 


(C.11) were the result. Arranging the equations as 


Wr - ce = antn + bn (eZ 20) 
dh 
Be att SLPO = a Ree) if Vea ce 
+ V1 (H,_j7h,) (C.11) 


allows a matrix equation to be formulated. The equations to 


be solved (analogous to C.1l) for the special stages were as 


follows: 
feed tray 
dh, 
WT. Sees QLP, = L,(h,-h.) os V, (h, -HL) 
+ Nese UH Meet ) + P(h,-h,) 
reboiler 
dh 


No enthalpy balance was required for the condenser. 


The matrix equation formed was 
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I< 


2 


As can be seen, terms in liquid and vapor flows were gathered 


on the right hand side, others on the left hand side. 


To solve for the liquid and vapor rates, the Matrix "A" 


is inverted. 
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Bbotn the MA" ematrix’inversionand the transformation of 
the "F" vector by the inverted A matrix were accomplished by 
computer library subroutines (IBM routine SIMQ). The result 
of these operations is the flow vector, which consists of the 


liquid and vapor flows for each stage. 


The second step in the calculation cycle is the integra- 
tion step. The composition balances were worked into the 


mers corm (C.14) 


xXx = 


\\00 


x Reapers ky 


To obtain this form, the following expression was sub- 


stituted for the vapor compositions. 


: (C.15) 
Yn ee 


The ratio KO was calculated for each stage just prior to the 
integration step, solely for the purpose of converting to 
the matrix form (C.14). As in the energy balances, the total 
balances were used to substitute flow rates for the holdup 


derivative term. The resulting set of equations have the 


general form (C.16) 
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dx 1 
dt wr (eins ae _ RX n * Os Se Soot 
dx 
nce ae) + xe (V_ (eK en ee 
at wre na, ond n n n+l n= 1 
+ i he See (G2L6) 
From Equation (C.16), the tri diagonal matrix B of order 


wr wiere "n" is the number of actual stages, can be repre- 


sented as: 


(Cork) 


IO 
i 


The parameters are given by: 


= = = - n = 2,¢ and<6, 9 (C228) 
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a, = wr (Y (1-K_) iD —) ve) 
a = aba. V (1-K ) 2 im) 
n Wr on n a 
b = iy K ) n=2,10 (Czae9)) 
n WT 1a HEY ees f J 
Coss GE ) n=199 (en20) 
n WT n+l d ‘ 


Upon inspection, this matrix can be seen to be diagonally 
dominant, as Ka the ratio of vapor composition to liquid 
Eomoosition on any tray, is always larger than 1, and thus 


oulhae a, terms are negative, and relatively large. 


To integrate Equation (C.14), the following sequence 
was used. An approximation to the derivative at time step 


K+l may be written 


cote es x ‘DT xk 


or x a oon ie x) /DP (eas 


where DT is the width of the time step. Substituting Equa- 


eron-(C.14) 
Bx aS < ais = xS) spr 
solving for cone 
Meble= (x.itpr Bayes xs (cee 4) 
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A better approximation for the derivative vector during 
the time step may be expressed as an average of the derivative 


at the beginning and end of the time step. 


: EGLO K+) eK 
xy = x + 5 x (Cr 289 


Meeecituting Equations. (C.14) and (C.21) 


ae S x) spr ag 5 span ae iE 5 aX xh 
Beaving for ai 
yktl = (I - pt/2 Be peat + DT/2 BX) x 


In this expression at time K, the matrix B at time K+l 
is unknown. If one assumes that the matrix does not change 
Suickly during the time step (i.e. if the system is not very 
non-linear), it may be taken to be constant. This results -in 


aL 
Bewation (C.24), from which xh 1 can be calculated 


eo Sy -ipt/2 Bs] a[P it r/2 Bax (C.24) 


Again, as in the coefficient determination step, library com- 


puter routines were used to do all inversions and matrix- 


vector operations. 


A simple Euler method was used to integrate the total 


mass balances, to obtain holdups. This set (3.30) of equd- 


tions was not able to be formulated into a matrix form similar 


to Equation (C.14), necessitating the use of the Euler or 


Other such method. 
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C.3 HOLDUP CORRELATIONS 


a) Rademaker et ite 


According to Rademaker at least two ways of expressing 


a holdup relationship can be found 


M=M {L, Up aca 
Or 


L 


L {M, Woes 


Rademaker presents parameters for these formulations. 
ments work, the Latter formulation was chosen, and a con- 


version was necessary. According to Rademaker: 


fe = eee) Mou [AV = x (40) ] 


as constant pressure was assumed, the pressure correction was 


dropped. Rearranging, 
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In the Rademaker formulation, 


6) 
au = ($=) = Zul Sec 
Vv 
§ 
ei 6s ($2) = 027 
M 


substituting these values 


AM, = PL IL TEING ie - aoa eV 
a i-l 


Since this is the perturbated equation, modifications 


must be applied for absolute values: 


M = M + MM, Ui = eA, VV) JS5 Ve OveoAV 
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z 
b) AICHE Bubble Tray Design Manual*” 


To predict holdups for use in the AICHE efficiency pre- 


diction procedure, the following equation has been determined 


empirically: 
iam Peo oet On Lows —0.05F + 0.0202 

where 

w = weir height (in) 

ioe = sli quid® flow’ (gpm/ft of flow width) 

i. 0<5 

F = Vola 

Vo = vapor rate (ft/sec) 

Og = gas density yeas 

eli height of clear liquid on theytray. (27) 


c¢) Francis Weir Formula 


; 36 2 
The expression given by Luyben™ 1s 
Paes. 332 (n) 3/? 
L 
ak, 3} 
where Bran I@quid rate (ft /sec) 
2 = length of weir (ft) 
h = height of liquid over weir (ft) 


while Perry gives the following expression 
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q = 0.415(L - 0.2 ho)h 
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A, = head above weir (ft) 
q = flow over weir (£t7/sec) 
g = force of acceleration due to gravity (£t/sec*) 


If it is assumed that 0.2h_ << L i.e. small head on weir, 


then the expression simplifies to 


feo 4ist ho /2g 


Peererncroducing that g = 32.2 ft/sec’, reduces the equation 
Bomene one cited by Luyben. Expressed in SI units, the equa- 


tion becomes 


ape Pee vias 


a{cm>/sec) 35 ..,4.. 5.8 (Cm) ho (cm) (2g(cm/sec )) 


Holdup values calculated at the standard operating conditions 


are given in Table C.1l for the Francis Weir and AICHE corre- 


lations. 
TABLE FCs 
CALCULATED TRAY LIQUID HOLDUPS 

tra Francis Weir AICHE 
Holdup (kg) Holdup (kg) 
it DELS iM} 
2 PLS 1.85 
3 2207 Poa 
+ 1.99 ‘ hemes fe 
5 1.92 1.65 
6 ioe 1553 
7 1.82 LoL 
8 L279 TAG 
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These holdups were calculated for the base operating state as 


follows: 
feed flow 18 g/s 
steam flow 14 g/s 
reflux Lilow 14 g/s 


Simulations using these operating conditions gave liquid and 
vapor flows for each tray, which were then used along with 


fray dimensions to calculate the entries in Table C.l. 


It may be noted that holdups predicted by correlations 
are much larger in magnitude than those measured. The AICHE 
correlation was developed for use with industrial scale 
columns so the pilot plant tower may not fit perfectly in the 
range in which the correlation applies. However, an indica- 
tion of holdup dependence upon column parameters was obtained, 
showing weak dependence upon vapor flow in one case and a 


weak dependence upon liquid flow in the other. 


C.4 EXPERIMENTAL PARAMETER SAMPLE CALCULATIONS 
eA] AICHE’ foldup Correlation 


The relationship for liquid tray holdup given by the 


AICHE Bubble Tray Design manual is: 


oat 1.65 + .019w - .65F + .020L 
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where F = fast te 
G 
Ze Netqne sor eligquidson the tray (in) 


V = vapor flow (ft/sec) 
QO. = vapor density (lbm/£t?) 
P= liguidetlow (gpm/ft of flow width) 


w = weir height 


mesumotion —- feet of flow width = 8 in = .7 £t 


example stage 3 (tray 2) 


oe 


Meee oos (liquid composition) 
y = .680 (vapor composition) 
v(£t>?/lbm) = = ize t 
for MeOH 
Beccdoex)) so82l acm 8, | ee | (454 _gmoles, (_lbmole 
fee) Gl atm)“ gmole °K © ‘28.322 Ibmole 32.04 Lbm 
w= 14455 ft? /lbm 
Similarly for water 
v = 26.83 £t>/1bm 
at y = .680 
v = .680(14.55) + .32(26.83) 
= 18.48 £t>/1bm 
gas volumetric flow is: 
bm) (18.48 ae 
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ual): one oa a Peete t racial Oa 
> ewe) volt. 
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area for gas flow: 


2 Z 
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Ean” ete 
aie ft/sec 
442 £t° 
= 1a 5. 9), st Ee, SeC 
gas density: 
-i 
V 
a ae ae 3 
= ee: iein/ fe. = 2054) lbm/ rt 
0 = (.0541) 
vapor term in correlation: 
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ee > 2139) 0541) °° = 324 
The liquid flow term: 
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cm ON Piss 3600 sec 
pee co co g/ SCC (33 gms’ (39-5 = athe ) 
mess its Lt 
= ,5625 gpm/ft of flow width 
ae Meo el ow sO obee 0 20G 


Meooete ado 42) = "65 0.c24)iae (20 ts5625) 


ee 
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but the column has 64 in of area 
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Mead density at x = 4.30 is .93 gm/em* 


Total weight of tray holdup = 1774 gms 
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eea.e Francis Weir Correlation 


The Francis Weir correlation in SI units is; 


q (cm/sec) => SCY ey Cored ho"? (cm) V2g (ence 


assume the weir is 3 sections of 2.5 cm widths 
example stage 3 (tray 2) 


muta £Low - 23.2 g/sec or 7.73 g/sec through one section 


Of the weir. 


mraqugad density at. .xi= .335: 
3 
or = -o 3° 4/7 cm 


volumetric liquid flow: 


3 
Geog. 73 9/sec)7(-93 g/cm?) 
= 8.3 om>/sec 
Substituting 


Seo se. 45) (2.9) no? °v98T 
h S209 
h = .402 


5 cm + .402 cm 


Total head on tray 


5.402 cm 
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total tray area 64 in? or 413 om” 


total holdup eae ON, Cas. Cm) 9s gm/em*) 


2074 gms 


C.4.3 Foss and Gerster Holdup Gorreiatinn © 


The Foss and Gerster correlation is: 


0.7.3 


fea 02 45702725 nh. - .29h Vio eae) 
W w aG 7 


where ho outlet weir height (in) 


i = gas velocity (ft/sec) 
Pg = gas density (ipmy fe) 
Geer ticiid f£low wate (cu) ft/sec) 
Z = average liguid flow width (ft) 
hy = liquid head on tray (in) 
@ssume average liquid flow width = 8 in or = ./ ft. 


example stage 3 (tray 2) 


Vos 15.1 c/s 
i252 g/7/SeEC 
1d eS a) 

Vo = 0 90 


from the sample calculation of the AICHE holdup prediction: 


ee 


UaQhoe = 7224 


volumetric liquid flow: 


2 

1_cm _ ift 

Gase2sn2 g/sec (—33-gms° (39.5 cm 
= 8.79 x 107? ft*/sec 
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De eee F067 25(2). = 229 (2) (.324) 4+ 448 000879) /(. 7) 


ie 0 Ie Or 3. 83° cm 


total mass holdup: 
Pee smem (4b3\em (.93 gm/cm?) 


= 1470 gms 


29 
eee. oc) ALCHE Tray Efficiency Calculation 


The correlation to determine efficiency given by the 


AICHE Bubble Tray Design manual is: 


eS 
Pee SP LoL+ .. 116 zo BPIEKS, Vea acO er cee Ls 
G y0- 
sc 
where 
'G 
Nee 5 (dimensionless) 
: Nene 
ie a height of outlet weir (in) 
Ug = gas viscosity 
5 
Pg = gas density (lbm/ft’) 
ie eeelTould race (gallon/min)/ft of average liquid 
“ flow width 
Da = diffusivity of the gas 


~ 
— 


assume average liquid flow width = 8 incon ee 


Limitations 
This correlation was developed from empirical data in 


the following ranges: 
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example stage 3 tray 2 


eee. Gg /S 
Ve= 7 t>.l g/s 
x7>..335 
y = .680 


Calculate the Schmidt number: 


Mee 3 
Meteamate the diffusivity using correlation from Perry: 


3 
P Tio lp 


Going through the steps of estimating atomic Taqia Etter. 


diffusivity was calculated. 
2 3 Z 
pe w0417 in /sec’= .269 cm» /sec 


vapor density: 


using ideal gas law, 


3 
VMeoH Se MRS tarsi y Ber 


from steam tables, 


Vv = 26.83 ft°/lbm 
water 
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vapor viscosities: 
from Perry, viscosities are 


water - .0114 cp (gm/cm sec x oa 


MeOH - .0116 cp 


assume an average of .0115 cp over all compositions. 


ge: oUAO00Droeomsemesesc se lE50 ar 


Noo 7 - 2 
PGD, 7209 ecMu/Asec 
= .487 
wer? = .70 
sc 


Liguid flow rate: 


3 
lcm Le 3 ,3600sec, 7 247 gpm 
e252 s (SS) (ee) (eS) ( ) 
La g/S (~o3Gm) ‘30.5em Ar ye 
= .395 gpm 
or .564 gpm/ft of flow width 
Vapor flow rate: 
from AICHE holdup prediction calculation, 
oD ae 
Vea aoe 24 
To calculate vapor transfer units: 
gs) 
sfiOne. 116.245 0.290V.9, + 0.02171, 


N = $i 
G e 
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Lid 
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OMe a6 (2) i= 0290". 3:24) +7 0297 (2564) ) 7.70 
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Calculate number of liquid transfer units: 


PIGS Le ey bag Pee yc. 


Oh) 
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Zz 
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+ 1.5) t, 


where D Mequid@-ditfusivity (£t7 hr) 


L 
feo 37.4 Z 62, 
L L 
W 
Zh Fe holdup, inches, of clear liquid on? tray 
Zz, = liquid flow distance (ft) 


From AICHE holdup sample calculation, 
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he Aes 
Vea = .324 
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apt Owe tl OS) at 
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85.1 sec 


6 Keane cm*/sec (Perry) 
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Liquid phase diffusivity De 
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Calculating the liquid transfer unit, 


55) 0 pedene saa) + 91.5) 585.8 
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Combine the liquid and vapor phase 


meg Sy 
= mG 
where XK = L 
m = slope of equilibrium line 
G = vapor rate Weies/ cea 
L = liquid rate Lbmoles/ft“hr 


From equilibrium data, at stage 3 composition 


m= .7064 
f= 923.2 g/sec 
Vea lS. kig/sec 


Assume approximately similar molecular weights in both 


liquid and vapor phase, 
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To convert to Murphree tray efficiency: 


eavculate eddy diffusivity 


os = (0.0124 + .0171 Ve +2002 5015 ee .0150w)- 
where Vo = vapor velocity (ft/sec) 
fe laguid bow (gpm/ft of flow width) 
w = weir height (in) 


substituting in precalculated values, 


(0.0124 + 0.171(1.39) + .00250(.564) + .0150(2))° 


D 
= 


.00456 


Calculate the Peclet number 
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00456) (85.2) 
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From Figure 8a, AICHE Bubble Tray Design Manual, a ratio is 


obtained from the Peclet number: 
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or Murphree vapor efficiency = .746 or 74.6% 
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APPENDIX D 


THE PILOT PLANT DISTILLATION COLUMN 


The following is a detailed schematic of the pilot plant 
mistillation column, with all pertinent control and data 


accumulation loops marked. 


Detailed information about all data accumulation or con- 
Heol loops on the column can be found in the DDC loop re- 
@ercs in this Appendix. DDC loops not marked on Figure D.1, 
but included in the DDC loop list are history data accumula- 


e1Oon loops: 


Engineering drawings and dimensions of the pilot! plane 
column and bubble cap trays are obtainable from the thesis 
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FIGURE D.1 A Detailed Schematic of the Pilot Scale Distillatio Colum. 
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LEGEND 
- indicates measurements available to IBM 1800 
LPID - DDC loop identification code in hexidecimal 
DAC - crosshatching indicates no loop 
ADC - multiplexer input address (in hexidecimal) 


DAC - multiplexer output address (in hexidecimal) 
- crosshatching indicated no output 
loop is for data acquisition only 


- indicates local analog equipment 
TYPE - gives function 


LYPE CODES 


T - transmitter 
R - recorder 
Cii=s controller 
Leas, indrcator 
Foi =) it LOW 
A - concentration (analyzer) 
 T +-£ temperature 
P - pressure 
DP - differential pressure 
n - identification number 


HXn | - heater exchanger 
n - identification number 


- gas chromatograph 


(z0) - thermocouple 


- indicates utilities 
UT - ST - 60 psig steam 
- CW - cooling water 


- solenoid valve 
Q n - identification number 


—<—_—-— - control lines 


wees ~ process lines 
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APPENDIX E 


THE GAS CHROMATOGRAPH ANALYSIS SYSTEM 


In this work, the bottoms composition analysis apparatus 
was changed from a Beckman series 4 industrial chromatograph, 
in which a common oven volume was shared by the separation 
column and automatic slider sampling valve, to a Hewlett 
Beemard gas chromatograph (G.C.) Model No. 5225A , with a 
separate oven for the sampling valve. Several problems were 
encountered during this transition. The approach to solving 


these difficulties is described herein. 


A short description of the steps involved in an analysis 
iS in order before discussing some of the problems encountered. 
A representative sample of the liquid to be analyzed must be 
collected and injected into the column. Since the bottoms 
product is a liquid, much smaller volumetric samples need be 
drawn from the process stream, in Vthe <order -CLPL microl 1 ere, 
which is about one tenth of the volume that would be used 


for a vapor sample. This creates two problems. It is more 


difficult to consistently provide a sample of the same size 


and furthermore a consistent analysis for small sample volumes 


becomes much more difficult. As a consequence of the small 


volume, the valve is prone to plugging by any particulate 


matter, and thorough filtering is necessary. 


i i : et 
- Sed 
a I ie i 
, ~~ = wu oo \ 
t . ! 
6° db i" -= 
Si i) 
>  WCaA fgiht , j 4 . ) & 
ath 2. ld _ 
.. ‘~~ ¥y 
<< i] 
ns : 
4 i : 
Pi. AN 
’ ~ F a ( 
~~ a a“ 
i 
i « 
—— 
» io - 
a } 


§ he -~ 
J 5 > 
> 
44 Peer & =~ 
#r- _ 
f ' y * 
= é “ 
- ae - 
hae 7 
* » ¢ @& 
44 a 
Se iu 
J - 
- 
- 
: = * a 
J | f 4 . ‘ Le | = 
4 - 
S 
- ¥ e 
- 
r= , z 
< oe, Bez Med apr Ps =» - ey 
y 44 pa4 © — alles 
7 ' 7 } 7 4 
e ie a? rh 
ie ' 
t » i . *\ > t 
- te t 
é a” r PoP \ yl 7» f 4 


: Pees t re : ; ‘ 
i! ] tg  panefop ‘edd awaoi.Be roe(el i 
. . | Pee hea th ae er -K: 
re mek Met Ppt! ace Pe Sa 
; 5 . ; a0 ¥ Ae tL J » € ban. : 
7 y = 


187 


A liquid sample must be quickly flashed before being 
introduced into the separation column. This required the 
sampling valve to be operated at higher temperatures than 
vapor sampling valves, and resulted in dissolved solids being 
deposited in the sample volume, again plugging it up. All 
these problems were encountered at some time during this work 
and three different sampling valves were tried and tested 


before the problems were sufficiently resolved. 


Pracerecne Most difficult of the above; problems to over— 
come was that of high temperature and proper flashing. Sche- 
matic diagrams of the valves used are shown in Figurs E.1l to 
B.4. With the arrival of the new gas chromatograph, the Valco 
valve shown in Figure E.2 was installed. Flow through the 
valve was poor, introducing a further time delay in the pro- 
cess measurement, which was undesirable. In addition, the 
oven in which the valve was operated initially had the Hewlett 
Packard insulation only, which was relatively poor, and slow 
flashing was the result. High temperature operation caused 


the valve internals to scorch, making actuation difficult. 


The internals were once changed, with little improvement. 


Other problems also surfaced but plugging, low flow and high 


temperature scorching were the most serious. A 10 micron 


filter installed in the sample line eliminated much of the 


plugging, but no amount of G.C. and valve oven temperature 


adjustment could avoid the scorching problem, and maintain 
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FICURE E.1 Beckman Liquid Sampling Slider Valve 
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good analyses. 


Another valve installed as shown in Figure E.3 
and E.4, manufactured by Applied Automation and known as 
Mark IX, contains very few moving parts and has a high metal 
content. This resulted in good flashing and negligible mecha- 
nical failure. The only problem encountered was a deteriora- 
tion of a teflon membrane with continued use. Adequate docu- 
mentation was not supplied with the valve, so after a trial 
and error procedure, proper operating conditions were found. 
A high temperature membrane was required, as the valve was 
run at 150°C, a temperature that was a compromise between 
deterioration and quick flashing. Valve actuation pressure was 
an important variable, minimal pressure being optimum. The 
internal arrangement of the valve, containing springs to load 
the flow paths was found to be important. A spring arrange- 
ment resulting in minimal pressure on the diaphragm is meces= 
sary for diaphragm longetivity. Other spring settings may be 


used for higher pressure, lower temperature operations. 


The internal flow path is also critical to good flashing, 


and good separation. The valve allows for a choice between 


two possibilities of the flow path. The one which results 


in better flashing was established by trial and error. A 


scheme was found to pass the liquid sample from the process 


stream to the separation column as a4 vapor, with consistency, 


reliability and reproducibility necessary for continuous 
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accurate analyses. It was then the task of the gas chromato- 
graph to separate and analyse the sample, and the task of the 
computer to interpret the G.C. output and calculate a compo- 


Sition in weight percent methanol. 


Problems encountered in the separation were typical of 
any gas chromatography analyses: the specification of ope- 
rating conditions for optimum separation and minimum elution 
times. Poropak Q packing had been used to separate water 
and methanol, but at small concentrations, the methanol peak 
was always found high on a long tail from the water peak, 
and made analyses difficult. Experiments were done with 
Poropak QS, a new packing, but it was found to degrade very 
Bapidly with time. Finally, the operating conditions listed 
in Table E.1 were established, based on a compromise between 


elution time and separation. 


TABLE Bl 


G.C. OPERATING CONDITIONS 


Separation Oven Temperature 30°C 
Valve Oven Temperature L502,0 
Carrier Gas (Helium) Flow 0.9 ml/sec 
Detector Temperature 200°C 
Detector Current 150 mA 
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A sample size of 2 microlitres was chosen, since for 
small concentrations of methanol, use of 1 microlitre sample 
volume gave rise to a methanol peak so small that it was lost 
in the noise encountered in transmitting the G.C. output 
signal to the IBM 1800 while use of larger samples tended to 
flood the separation column and degrade the separation. It 
should be noted also that subsequently an amplifier was in- 
stalled on the G.C. output, thus reducing the significance 


of noise. 


The methanol content in the samples was determined by 
the G.C. analysis package on the IBM 1800. This program 
integrates the area under component peaks, and by comparing 
relative areas, relative quantities in the sample are deter- 
mined. Operation of the program is by specification of a 
"GC job" for each particular type of separation or chromato- 
gram. Peaks beginnings and ends are identified by observing 


first and second derivative changes in the G.Cc. output. De- 


fining an appropriate G.C. job for a typical chromatogram is 


accomplished with the help of a new GCRAW group of programs 


which collect and plot peak first and second derivative data. 


The addition of the amplifier was the single most impor- 


tant modification that reduced the di fticuleyyin obtaining 


consistent measurements, once the sample had been injected 


into the separation column. The significant changes in the 


signals have already been discussed in Chapter 5. 
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Degradation of analyses can be detected by a change in 
the elution times of the peaks. Once the G.C. has warmed up, 
elution times should not vary by + 2 seconds. Flow through 
the sampling valve is a reliable indicator of valve failure. 
Inconsistent analyses may be traced to valve failure in most 
cases. The line filter must be cleaned periodically and care 
must be taken to keep the valve oven temperature as close to 
the set operating temperature as possible. Also, since the 
column is operated as a closed system, any solids accumulate 
in the methanol-water solution; the solution must be purged 
from time to time, or plugging of the valve by dissolved 


solids will become a problem. 
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APPENDIX F 


EXPERIMENTAL PARAMETER PREDICTION 


The parameters of holdup and efficiency are important 
in determining the agreement between experimental and simu- 
lated distillation column responses. Experimental measure- 
ments were performed to estimate parameter values, or to 
define a parameter range. The results of these measurements 
are documented in this appendix. In addition, parameter 
values calculated from several common correlations are also 


presented. 


Pyl Sample Calculations for Experimental Holdup Measurement 


First estimates were made of relative volumes of liquid 
on all trays, using tray 5 as the reference’ tray. PFaverin— 
dependent estimates were made, with the column operating at 


steady state and the results of the different observers ave- 


raged. 


Sample (for illustration only, not actual data): 


TRAY u Z 


RATIO aoe Epes anak Heal un 
The column was shut down abruptly; all input and output streams 


were stopped. All liquid was drained into a container and 
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total weight of solution in column Z20.k6 
average composition (measured) 25% MeOH 
density ~952 KO/x 


volume DOK / oD 2) KO/ © 


= 2 08 


The measured volumes of the reboiler and condenser, being 
constant were subtracted: 
measured reboiler volume PcG 


measured condenser volume LW 


Then the total tray volume was calculated 


total tray volume 21.0 - 7.3 - Le TAS ae 


The estimated ratios were converted into fractwonsrand 
the total tray volume multiplied by the fraction to obtain 
individual tray volumes. An approximate density was applied 
to each tray, yielding the tray mass holdup, as given below 


in a sample set of data (not representing any actual measure- 


ments) : 
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Tray Ratio: Fraction Volume (2) Density Mass (kg) 
(kg/%) 
il lear ~ £45 ath 0 Aas FES 
2 ee 45 Wy E3 0 Leeks 
3 ee wis 9 ses) deo: 
4 vs, i RE ChS) 159 295 aces 
5 Abe e) 2.0 aa 220 nly. WE 
6 0.9 Woks pene 90 in Se 
I Or o> 114 Pes Backs iG 
8 OS) ya 02 ees «8D 1.04 


Total SAe3 
Eraction 12/3. 43-'= 2145 orstirsthieray 
A summary of all estimated ratios and the corresponding tray 


Meidups for a variety of operating conditions is given in Table F.1. 


fee, Attempts At Correlation of Measured Holdup 


With the nonlinear model, the steady state terminal com- 
positions were matched with their experimental counterparts, 
by fitting efficiency and heat loss parameters in the model. 
The model yielded also intermediate compositions and flow- 


rates, which were used as independent variables in a corre- 
lation technique. 


Simulated steady state values of liquid flow rates, 


vapor flow rates and corresponding measured tray mass holdups 


were used as input data to a correlation program on an IBM 


370 at the University of Alberta. The name of the 


Operating 


TABLE Fol 


EXPERIMENTAL HOLDUP DATA 


Conditions 
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correlation (regression) program was MLREGR and the desired 
format of the correlation is given below, where P denotes 


fitting parameter. 


The regression gave no meaningful results, due to the 
excessive scatter of the data. This data is presented in 
Figure 5.10, where holdup is plotted against liquid flow rate, 
and where it may be observed that there is no discernible 


pattern. 


It would have been expected that even using only liquid 
and vapor flow as independent variables that some valid cor- 
relation would have resulted. Perhaps the visual observation 
technique led to such variation that any trends were obscured. 
Even when tray holdups are plotted against feed flow rate 
suchi as in Figures F.1 to F.4, considerable scatter is ob- 


served, although a general trend is discernible. An increase 


in feed rate appears to bring about an increase in holdup on 


all trays. The holdup values in these figures represent a 


graphical summary of all values that could be calculated from 


the different independent observations. 


F.3 Efficiencies Prediction 


An experimental and a mathematical technique was used to 


Beeeablish a range of efficiency that would be used in the 
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TRAY HOLDUP (KG) 
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FEED FLOW (G/S) 


FIGURE F.1 Experimental Holdup Estimates, 
Trays l & 2 
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FIGURE F.2 Experimental Holdup Estimates, 
Trays 3 & 4 
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TRAY HOLDUP (KG) 


TRAY HOLDUP (KG) 


16 18 70 59 
FEED FLOW (G/S) 


Tray 6 


16 18 20 Fig 
FEED FLOW (G/S) 


FIGURE F.3 Experimental Holdup Estimates, 
Trays 5 & 6 
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FEED FLOW (G/S) 


Tray 8 


16 18 20 Ze 
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FIGURE F.4 Experimental Holdup Estimates, 
Trays 7 & 8 
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Simulation. 


At steady state, liquid samples were taken from each 
stage of the distillation column. The samples were subse- 
quently analyzed by a gas chromatograph and a column compo- 
sition profile was determined. This profile was matched by 
simulation with the mathematical model. With the same ope- 
rating conditions, efficiencies were changed by trial and 
error, to give new steady states. The efficiencies that gave 
an initial steady state composition profile within a few per- 


cent of that measured are shown in Table F.2. 


This method has two disadvantages. First, the accuracy 
of the mathematical model must be accepted, before actual 
comparison to experimental results are carried out. Secondly, 
it is a time consuming and expensive procedure, in terms of 
computation, since the values must be established by trial 
and error. If vapor compositions could be sampled and mea- 


sured, efficiencies could be calculated directly. 


The only prediction attempted was done using the recom- 


mended AICHE efficiency prediction for bubble cap trays. 


This method has limitations (see sample calculations, Appendix 


C) for pilot plant scale operations, and the calculated esti- 


mates are included in Table e.2.° The Peciet numbers offer an 


Hadi cation of the extent of mixing that tale place; a Pe of 


near 1 indicates good mixing, while a higher Pe of greater 
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than 5 indicates an approach to plug flow. The buble cap trays 


here all appear to be well mixed. 


mea. Estimation of Column Heat Loss 


As presented in Appendix B, the programs DASS, DATAL 
and BALNC calculate an energy balance on the dvsti it auLon 
column, with the difference between energy input and removal 
taken as heat loss. This heat loss was found to be due to 
losses from the column walls and losses from the gas chroma- 
tograph recycle system. The combined heat loss as calculated 
by the DASS program was of the order of 8-9 kJ/s. To verify 
these calculations to some degree, the following heat loss 


calculation was performed. 


a) column 


7 trays - 23 cm diam x 30 cm height 
l tray - 20 cm x 60 om + 9 cm x 120 °cm 
lL reboiler + glass below first tray = area of 


11/2 regular trays. 


heat loss due to convection (45) 


. > : t 
q=h Bit-T 9 where ee is the ambient temperature 


assuming natural convection, laminar flow around column, 
: 2 
for vertical cylinders 
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Using the temperatures predicted by the computer simu- 
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the nature of the correlations, 
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heat transfer coefficients are calculated. 


Due to 


calculations were done in 


Engineering units, and results were converted incor Si Wace a. 
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TABLE F.4 


HEAT LOSS DUE TO RADIATION 


Stage 7g uy Oaare 2 2 PTOTAL 
(°R) (°R4) (°R4) (feo (joules/sec) (Kj/s) 
(x1011)  (x1011) (x 1077) 

dl 672 2.04 Tae PA 5 Oe BZ0 
u 665 e396 ei O7) ey, 050939 omy 
3 656 1 CO LOG dea Fe 02085 PERS 
4 647 Lees 0.96 aye) O07 7 a 
S) 638 E266 O67 Let O60 eke 
6 oleh Maes) 0.80 LT ye 0.064 Swe 
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8 616 1.44 Oo ey), O05 09 
9 609 The Bes: Oe) 4.90 Opa E STS, Ree: 


Note that the total of convective and radiation heat loss 


(the sum of the last column in Table F.4) explains one quarter 


of the steady state heat loss. 


Bb) G.C. recycle system 


Heat loss was not rigorously calculated here due to the 


complex configuration of the recycle system; the circuit 


included a small heat exchanger, 4 pump, lengths of tubing 


of different sizes and materials and various other points of 


Mee raet with metal surfaces in the system. Instead, the 


material and energy balance programs DASS, DATAC and BALNC 


Mee en with and without the recycle system operative. From 
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mrese runs, it was found that half of the total heat loss 


could be attributed to the G.C. sampling system. 


Meo) energy Storage in Column Material 


As pointed out by Rademaker et Buh the storage of 


energy in the metal and glass of the column may be an impor- 
tant enthalpy consideration, more important perhaps than the 
energy storage in the liquid holdup on the tray. Tis: ‘ettect 


was introduced into the model in the enthalpy balance. 


Volume of metal in a tray: 


plate itself 3/38" x 14" diam (less cuc- 55.00 (osmouy 
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downcomer Ve aOR a nny eo ke Oe 234 
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TOTAL 655 51ne 
0378 a 
Total volume , 4 
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Stainless steel density 


Total weight 18.45 lbm 


° 
Heat capacity 0.11 BTU/lbm°F 


To raise tray metal one degree Celsius, the amount of heat 
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mimelariy, Lor the tray glass, 


Setamex L2') x, 1/8" SRA ine 
total volume 0218 £t> 
glass density 170 ien/ ee 
heat capacity 0.2 BLU ALom? Fr 


Energy storage per °C is: 


Q = mC 
2 Pp 


A peibm. (0.2 BIU/lbm°er) (1.054 Rj/BTU) (2.3) C7) 


1.58 Kj 


Energy storage per tray and glass secticoneper “°C is thus: 
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APPENDIX G 


EXPERIMENTAL STEADY STATES 


The pages following present a summary of the terminal 


operating states of the distillation column open loop 


responses presented in this thesis; the run number codes 


are explained below: 
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SrecaoY STATE CONDITIONS BASED ON 
RUN NOQ OLIOSE 


PEED FLOW 
pee LUX, FLOW 


SEAM 


BOTTOM PROD 
TOP PROD 


COOL WATER 


FLOWW = 


DIST COMP = 
BOTTOMS COMP= 
FEED COMP = 


PRESSURE 
COND LEVEL 


REBIR 


DiFe PRESS 


REBIR 
PLATE 
PLATE 
PEATE 
PLATE 
PLATE 
PLATE 
PLATE 
Peg a 


COND TEMP 


STEAM 


COND* T TEMP 
mee Lux eb EMP 
Reco PeEMP 

BOTTOMS TEMP 


REB'R 


BeeD INLET 
BEeLUX INLET 
eat ON HEAD 


WATER 
WATER 


BEVEL 


O'HEAD= 


aNown FWN re 
o 
m 
= 
0 
" 


TEMP 


mEMP 


INLET 
OUTLET= 


20 POINTS 
02/06/76 - : 
G/SEC DEV= 0.0368 
G/ SEC DEV= 0.1511 
G/SEC DEV= 0.2846 
G/ SEC DEV= 0.4370 
G/SEC DEV= 0.0895 
G/ SEC DEV=18.0728 
WTZ MEOH OEV= 0.0791 
WT% MEOH DEV= 0.0000 
WT% MEOH  DEV= 0.0000 
KPA DEV= 0.6574 
CM DEV= 0.0410 
CM DEV= 1.7406 
KPA DEV= 0.1655 
DEG C DEV= 0.1393 
DEG C DEV= 0.1238 
DEG C DEV= 0.1237 
DEG C DEV= 0.1295 
DEG C DEV= 0.1295 
DEG C DEV= 0.1189 
DEG C DEV= 0.1168 
DEG C DEV= 0.1026 
DEG C DEV= 0.1080 
DEG C DEV= 0.0993 
DEG C NEV= 0.1333 
DEG C DEV= 0.1416 
DEG C DEV= 0.1351 
DEG C DEV= 0.1023 
DEG C DEV= 0.0869 
DEG C DEV= 0.1439 
DEG C DEV= 0.1617 
DEG C DEV= 0.3211 
DEG C DEV= 0.1026 
DEG C DEV= 0.0895 
DEG C DEV= 0.2091 
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BOTTOM PRODUCT 
TOP PRODUCT 
ESOSURE ERROR=PC 


ShEADY «STATE. 2 
RUN NO OL106B 


05/06/76 
Pos3s6 ASSES B80T 
ia 22 9G4oec ee 
L4696G/SEC COO 
oe Bo 
96<50WT4 iMEOR BOT 
LeeoDEG € Ree 
LO7«4DEG C PRE 


Mat T Gua yk eo BA 
FLOW 
LG7SEC) hat 
16.99 
9.06 
8.12 
el 
E NaGeR eGaY © 6B 4 
ENTHALPY 
[Ka/ SEC 
COOLING WATER 23-87 
REFLUX 2463 
TOP PRODUCT 
STEAM eke: 
BOTTOM PRODUCT 


HEAT LOSS 
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TOM PROD 920667 SEG 
PROD < aMule Cys ate: 
L WATER 287.464G/SEC 
DeGeMe 50 .O00WT% MENH 
TOMS COMP 3-80WT3 MENDH 
bu a DNOE 62.9065 08 
SSURE ~28-9KPA 
Lasley Ni LGus 
COMP METHANOL WATER 
WieukCiWwis (GGASEC) (G/SEC) 
50.00 8-49 8.49 
3-80 0.34 Set 
96.50 he 83 0-28 
-—3.6 ae, 
[aaa NC rs 
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y (SS eG) 
47.59 
mars gi 
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3,55 
61.88 
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STEADY STATE CONDITIONS BASED ON 
RUN NO OL106B 


FEED FLOW = 16.990 
REFLUX FLOW = 13.596 
STEAM FLOWW = 14.967 
BOTTOM PROD = 9.062 
TOP PROD Sig kG elec) 
COOL WATER =287.641 
DIST COMP = 96.509 
BOTTOMS COMP= 3.800 
FEED COMP = §0.000 
PRESSURE =-28 .931 
COND LEVEL = 17.166 
REBtR LEVEL = 40.636 
DIFF PRESS = -5.913 
REB'*2 O'HEAD= 93.9 
PLATE L TEMP= 83.4 
PLATE 2 TEMP= Wil ee 
PLATE 3 TEMP= 114.0 
PLATE 4 TEMP= 7332 
PLATE 5 TEMP= 68.5 
PLATE 6 TEMP= 66.0 
PLATE 7 TEMP= 64.5 
PLATE 8 TEMP= 63-3 
COND TEMP = 61.6 
STEAM TEMP = 10724 
COND'T TEMP = 106.3 
REFLUX TEMP = Spl 
FEED TEMP = 36.4 
BOTTOMS TEMP= 4761 
REB'R TEMP = 93.9 
FEED INLET = (Zin6 
REFLUX INLET= 63.9 
COL O'HEAD = Se. 


WATER INLET 
WATER OUTLET= 


19.8 
39.2 


50 POINTS 
05/06/76 
G/SEC NEV= 0.0473 
G/SEC DEV= 0.1908 
G/ SEC DEV= 0.1997 
G/SEC NDEV= 0.0583 
G/SEC DEV= 0.1266 
G/SEC NEV=14.4864 
WT% MEQH NEV= 0.1507 
WT% MEOH NEV= 0.0000 
WTS MEOH NEV= 0.0000 
KPA NEV= 0.6172 
CM NEV= 0.0888 
cM DEV= 0.2578 
KPA DEV= 0.2248 
DEG C DEV= 0.1889 
DEG C NEV= 0.1621 
DEG C NDEV= 0.1366 
DEG C NEV= 0.2070 
DEG C NEV= 0.1469 
DEG C NEV= 0.1456 
DEG C NDEV= 0.1403 
DEG C NEV= 0.1246 
DEG C NEV= 0.1238 
DEG C -DEV= 0.1374 
DEG C NEV= 0.2012 
DEG C NEV= 0.2012 
DEG C NEV= 0.1289 
DEG C NEV= 0.1035 
DEG C NEV= 0.1141 
DEG C NEV= 0.1989 
DEG C NEv= 0.1581 
DEG C DEV= 0.2741 
DEG C DEV= 0.1369 
DEG C DEV= 0.1111 
DEG C NEV= 0.2420 
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FEED FLOW 
REFLUX FLOW 
STEAM FLOWW 
FEED PLATE 
Beis. COMP 
seco INLET 
SecA) TEMP 


FEED 


BOTTOM PRODUCT 


TOP PRODUCT 


CLOSURE ERROR-PC 


STEADY STATE DATA 
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10.49G/SEC 
6. 1/1G/ SEC 
327.214G/SEC 


RUN NO OLLO6E 
05706776 
16-99G/SEC BOTTOM PROD 
LoOwlocsoce TOP PROD 
14.99G/SEC COOL WATER 
4 FEED COMP 50.-00WT% MEOH 
97.09WT% MEOH BOTTOMS COMP 17 -O0W le MERE 
TZ2.00EG. GC REFLUX INLET 63.4DEG C 
Ole 6 GeC PRESSURE -28.3KPA 
MalAl 7 o& Ry Tacks il 8 Avge A oN eee 
FLOW COMP METHANOL WATER 
(GASEC } (Wie PGi) uo soecy LGA SECH 
16599 showed 49) 849 8649 
10.49 166-00 Leo 8-70 
ere it 97 09 Sood Die ies 
lea -2.2 4.8 
NE yRaGe Y Biot Nass mean ome 


ENTHALPY IN 


KS/SEGE 
COOLING WATER 26.09 
REFLUX 3.09 
TOP PRODUCT 
FEED sae 
STEAM 40.26 
BOTTOM PRODUCT 

73.90 


TOTAL 


HEAT LOSS 


ENTHALPY OUT 
CXS ece 


Stirs iad 
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3.64 
65.8% 
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STEADY STATE CONDITIONS BASED ON 50 POINTS 
RUN NO OLLOG6E 037 0677715 


FEED FLOW = 16.990 G/SEC NEV= 0.0533 
REFLUX FLOW = 16.186 G/SEC NEV= 0.1587 
STEAM FLOWW = 14.993 G/SEC NEV= 0.1928 
BOTTOM PROD = 10.492 G/SEC NEV= 0.0899 
TOP PROD = 186007 el -G7SEC DEV = Oro Zi 
COOL WATER 2=327.146 G/SEC NEV=12.3664 
NIST comp = 97.099 WT% MEOH NEV= 0.1334 
BOTTOMS COMP= 17.000 WT% MEOH NEV= 0.0000 
FEED COMP = 50.000 WTY MEOH NEV= 0.0000 
PRESSURE =-28.878 KPA HEV=" "0.4081 
SOND “PEVEL = 17.283 CM DEVe “0.0521 
REB'R LEVEL = 43.577 CM NEV= O.4103 
DIFF PRESS = -4.372 KPA NEV= 0.1229 
REB'R O'HEAD= 85.8 DEG C NEV= 0.1758 
PLATE 1 TEMP= 15.0 DEG NEV=Oeleare 
PLATE 2 TEMP= TZ. DEG UC DEV =r aa) 
PLATE 3 TEMP= 116.0 DEG C NEV= 0.1967 
PLATE 4 TEMP= 71 sO DEG *€ DEV=70.1204 
PLATE 5 TEMP= 66.7 DEG C NEV= 0.1393 
PEATE 6 TEMP= 64.4 DEG C NEV= 0.1089 
PLATE 7 TEMP= ° 63.2 DEG C NEV= 0.1293 
PLATE 8 TEMP= S22 eee G NEV= 0.1088 
COND TEMP = 6150 "DEGAS NEV= 0.1091 
SUEAMe TEMP = 101-2 DEG C NEV= 0.1889 
COND'T TEMP = 98.7 DEG C NEV= 0.1907 
REFLUX TEMP = 51.4 DEG C NEV= 0.1250 
FEED TEMP a 36.7 DEG C NEV= 0.1063 
BOTTOMS TEMP= 47.4 DEG C NEV= 0.1093 
REB'R TEMP = as.7 Dee -C NEV= 0.1957 
FEED INLET = 72.0 ESAe FEV= 0.0502 
REFLUX INLET= 63.4 DEG C DEV= 0.3800 
COL O'HEAD = 62.4 DEG C NEV= 0.1189 
WATER INLET = 19.0 DEG C DEV= 0.0980 
WATER OUTLET= 37 <8 DEG Ss REV= Gacue? 
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fee) FLOW 

pee LUX FLOW 
STEAM FLOWW 
Beet) PLATE 
mest COMP 

eee TNLET 
aon TEMP 


REED 


BOTTOM PRODUCT 
TOP PRODUCT 
moe suURE ERROR-PC 


SHADY 5. VAbe eDAgA 
RUN NO OL1078 
067-05. /.10 


POwo) Gy Sree BOTTOM PRON 
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pe rg Cy na & 


13.47G/SEC TOP PROD 8.06G/SEC 
14.94G/SEC COOL WATER 398.11G/SEC 
a Febp= coup 50.00WT% MEOH 
96.51WT% MEOH BOTTOMS COMP 4.40NT% MEQH 
71.5DEG C REFLUX INLET So e7DECr 
105.6DEG C PRESSURE -30.2KPA 
Tepe eA LB AES cme 
FLOW COMP METHANOL WATER 
GG7SEC):: (Nae Grete Sec) (G/SEC) 
16.97 50.00 8.48 3.48 
9.12 4.40 0.40 Bi Fal 
8.06 96-51 7.78 0.28 
Ive 2 -3.5 6.0 
ete args YA. 8 Ae eee 


ENTHALPY IN 


(KJ/SEC) 
COOLING WATER 30.32 
REFLUX 209% 
TOP PRODUCT 
See) 4.4) 
STEAM 40.24 
BOTTOM PRODUCT 
ia eek 


TOTAL 


HEAT LOSS 


ENTHALPY OUT 
(KI7 SEC} 


eis wee ke 
2043 
1.45 


6.58 
Z< SU 
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STEADY STATE CONDITIONS- BASED ON 
RUN NO OLLO7B 


PEED FLOW 
REFLUX FLOW 
STEAM FLOWW 
BOTTOM PROD 
TOP PROD 
WATER 
es nC OMP 
SetTOMS COMP 
BeeD COMP 
aeos URE 
GONOD “LEVEL 
peo’ R LEVEL 
eer PRESS 


COOL 


REB'R 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE 
PLATE 


BOND <lEMP 
STEAM TEMP 
eee ND! Ta eT EMP 
eeriUXxX TEMP 
eee) TEMP 
BOTTOMS TEMP 
Ree’ R TEMP 
Seep INLET 
eer lUxX INLET 
Bok O'HEAD 
WATER INLET 


ONAN FWY rE 


YT Wig |e ceva 


un bu db db yp wy we bw 


WATER OUTLET= 


Mera chia 
13.476 
14.946 
Fm ou 
8-067 


Ge kel bie) 


26 eo kG 
4.400 
50-000 


ican > 
Bite o 
=5 +206 
OZ ek 
ot 
Tae) 
betes @) 
rae, 
67.2 
64.7 
6363 
2%. 1 
60.5 
105.6 
104.3 
S0R0 
35.0 
45.6 
Zoi 
(aes 
62%. ih 
S265 
Veec 
33.6 


50 
S67 OSG 

GASEC NEV= 
G/SEG NEV= 
GHSEG NEV= 
Gis ec NEV= 
GSEC NEV= 
G/SEE NEV= 
WT% MEOH DEV= 
WT% MEOH DEV= 
WT% MEOH NEV= 
KPA NEV=s 
CM NEV= 
CM NEV= 
KPA DEV= 
DEG iC DEV= 
DEG ¢ DEV= 
DEG. DEV= 
DEC RE NEV= 
DEGwG DEV= 
DEG BE NEV= 
DEG: /G DEV= 
DeGaes DEV= 
DEG ac NEV= 
DEG EG NEV= 
DEG AC NEV= 
DEG -¢ DEV= 
DEG: <C NEV= 
DEG DEV= 
DEG C NEV= 
Hie GaG NEV= 
BEG. Cc NEV= 
DEG-.G NEV= 
DEG eG NEV= 
DEG G DEV= 
HEG Cc NEV= 


POINTS 


0.0417 
0.1639 


0.2142 


0.0495 
Ore lG2 7 
8.6673 
0.1481 
0.0000 
0.0000 
0.3451 
0.0446 
O2525 
Ore FL 
0.1861 
0.1602 
0.1064 
0.1964 
0.1390 
On1s2o 
Oe E253 
0.1236 
0. 1236 
Grit 202 
Ore 20:12 
0.2204 
0.1471 
0.1010 
0.1070 
OSL699 
O.1520 
0.4650 
0.1242 
0.0873 
Oe LOuee 
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SEA go Pale aaa 
RUN NO OLILO7E 


05/06/76 

Beep PLOW £6 S9967S5EE BOTTOM PROD S2556/S51EC 
REFLUX FLOW POSEIG 5 Se TOP PROD S626G/75eC 
STEAM FLOWW L4.976/SEC COOL WATER 282-43G/SEC 
peep “PLATE 4 FEED SeONE 50.00WT% MEOH 
eS eocOMP 94-67WT% MEOH BOTTOMS COMP O.1L1O0WT% MEOH 
pee SENET GieoleGac REFEOeI NLS 62.9DEG C 
STEAM TEMP POO s3DEG C PRESSURE =30.2K PA. 


MGA hile O5A 2 Bo A ot YA" NiaGons 


FLOW COMP METHANOL WATER 

(G/SEC) (WEYECTL ING SEG) (G7SEG) 

Been [6.99 50's 00 8.49 8.49 

BOTTOM PRODUCT. Bie. 35 Biggie 0.00 8.34 

TOP PRODUCT Steud 94.67 S50 0.47 
BECSURE ERROR-PC 1.0 dines eh aii 


Eo Ne SREGEY 5: Ay DRAaN EG3e 


ENTHALPY IN ENTHALPY OUT 

Het SECs (KJ/SEC) 

COOLING WATER ei ep 43.95 

REFLUX 2206 1.99 

TOP PRODUCT 1.63 
FEED 4.42 

STEAM 40.41 6.87 

BOTTOM PRODUCT 3.40 

TOTAL ii ees eu 57.85 

Toe 


HEAT LOSS 


wit C1? 7 Ak 0, -— 


i ¥ . © ‘ a 7 
4 ; ' : 
~ ye) aN 
; } — 7 , va 
" iba nm 
is 
4 " Py a - 


, ar ac: oT hse yonere « 7a 
at 010 ‘eu uel <, 
BT \aONEH eet 


s 
t oa ; b 


~ ; 4 Pan 
a art oe 7: Ce Oli 
BORG DT Ja? TO) 


sarap Ja00 Ts NFP 
GwOA. 2MCVT'T IE iad vivties 


+ at WISE.» “F MGA, 1§ 
‘SRL Ry |, _ ‘ DIGEz OT 


rs, 
iv gid 08° a AO ae ch 


4° uM < MiyD — wn »= I p+ y J 
37iDb ae Pein’ (SREVES: . 


: y 
Pare, Alig OF 70 64 ene 
IO! Oe! J 6.9." Vip s. 
= ° 6.0" i £855 a aan. 


4 3 Ms ye is aia | wR ae. ws 
: ; % mane 
ruc ¥e jah Twa Al: as pe 
Lee Va bet 1G 1 pare aie ji . 
iba rae, ‘eran | i 
aarn eh ue 


eEeADY STATE CONDITIONS BASED ON 
RUN NO-OLLOTE 


FEED FLOW 
ReeLUA FLOW 
STEAM FLOWW 
BOTTOM PROD 
foe RROD 
COOL WATER 
en sa -COMP 
BOTTOMS COMP 
PEED COMP 
PRESSURE 
BOND. LEVEL 
wep rR LEVEL 
Doere PRESS 
ees R O' HEAD 
eae lt TEMP 


PLATE 2 TEMP 
PLATE 3 TEMP 
PLATE 4 TEMP 
PLATE 5 TEMP 
PLATE 6 TEMP 
PLATE 7 TEMP 
PLATE 8 TEMP 
COND TEMP 

STEAM TEMP 


COND? TEMP 
REFLUX TEMP 
e2eO sFeEMP 
BOTTOMS TEMP 
Reo TeMP 
Boece. PNCET 
REFLUX INLET 
COL O'HEAD 
WATER INLET 
Pare GUTLET 


nuuuw tb bon won wou wet 


nun uu ou bon bu 


50 POI iS 


NEV= 0.1384 
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DEG C NDEV= 0.1772 
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STEADY STATE CONDITIONS BASED ON 
RUN NO QLLOS8E 


REED FLOW 
REFLUX FLOW 
STEAM FLOWW 
BOTTOM PROD 
TOP PROD 
COOL WATER 
DS eC OMP 
BOTTOMS COMP 
SEE DEC OMP 
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GUND. LEVEL 
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REFLUX TEMP 
Pe een 
BOTTOMS TEMP 
REBeR TEMP 
PEED INLET 
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0.2495 
0.0600 
0.1090 
9.7596 
0.1545 
0.0000 
0.0000 
0.3768 
OLN OW ATE: 
0.3674 
Ora CON 
0.1646 
0.1684 
Geeore 
Geezer? 
0.1816 
epoe i Res) 
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0.1458 
0.1039 
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0.0874 
0, 1602 
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(G/SEC) (Wie eP.C Ty), sGGAS EC) V675 667) 
Peed ete 20 50.00 8-48 8.48 
BOTTOM PRONUCT 9.07 4.60 0.41 8-66 
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